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1. Introduction 

The U(1)a anomaly has played an important historical role in establishing QCD as the 
theory of the strong interactions. The description of radiative decays of the pseudoscalar 
mesons in the framework of a gauge theory requires the existence of the electromagnetic 
axial anomaly and determines the number of colours to be iV c = 3. The compatibility of 
the symmetries of QCD with the absence of a ninth light pseudoscalar meson - the so- 
called 1 U(1)a problem' - in turn depends on the contribution of the colour gauge fields to 
the anomaly. More recently, it has become clear how the anomaly-mediated link between 
quark dynamics and gluon topology (the non-perturbative dynamics of topologically non- 
trivial gluon configurations) is the key to understanding a range of phenomena in polarised 
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QCD phenomenology, most notably the 'proton spin' sum rule for the first moment of the 
structure function g^. 

In this paper, based on original research performed in a long-standing collaboration 
with Gabriele Veneziano, we review the role of the U(1)a anomaly in describing a wide 
variety of phenomena in QCD, ranging from the low-energy dynamics of the pseudoscalar 
mesons to sum rules in polarised deep-inelastic scattering. The aim is to show how these 
experiments reveal subtle aspects of quantum field theory, in particular topological gluon 
dynamics, which go beyond simple current algebra or parton model interpretations. 

We begin in section 2 with a brief review of the essential theoretical toolkit: anoma- 
lous chiral Ward identities, Zumino transforms, the renormalisation group, and the range 
of expansion schemes associated with large N c , notably the OZI approximation. Then, in 
section 3, we build on Veneziano's seminal 1979 paper [1] to describe how the pseudoscalar 
mesons saturate the Ward identities in a way compatible with both the renormalisation 
group and large- N c constraints and derive a generalisation of the famous Witten- Veneziano 
mass formula for the rj' which incorporates, but goes beyond, the original large- N c deriva- 
tion [2, 3]. 

In section 4, we turn to QCD phenomenology and describe how this intuition on 
the resolution of the U(1)a problem allows a quantitative description of low-energy pseu- 
doscalar meson physics, especially radiative decays, the determination of the pseudoscalar 
decay constants, and meson-nucleon couplings. We review the U{1)a extension of the 
Goldberger-Treiman formula first proposed by Veneziano [4] as the key to understanding 
the 'proton spin' problem and test an important hypothesis on the origin of OZI violations 
and their relation to the renormalisation group. Low-energy rj and rj' physics is currently 
an active experimental field and we explain the importance of an accurate determination 
of the couplings g^NN and g v >NN in elucidating the role of gluon topology in QCD. 

All of these low-energy phenomena have counterparts in high-energy, polarised deep- 
inelastic scattering. This enables us to formulate a new sum rule for the first moment of 
the polarised photon structure function g\ (section 6). The dependence of this sum rule 
on the invariant momentum of the off-shell target photon measures the form factors of the 
3-current AVV Green function and encodes a wealth of information about the realisation 
of chiral symmetry in QCD, while its asymptotic limit reflects both the electromagnetic 
and colour U(1)a anomalies. We show how this sum rule, which we first proposed in 1992 
[5, 6], may soon be tested if the forthcoming generation of high-luminosity e + e~ colliders, 
currently conceived as B factories, are run with polarised beams [7]. 

The most striking application of these ideas is, however, to the famous 'proton spin' 
problem, which originated with the observation of the violation of the Ellis-Jaffe sum rule 
for the first moment of the polarised proton structure function g^ by the EMC collaboration 
at CERN in 1988. This experiment, and its successors at SLAC, DESY (HERMES) and 
CERN (SMC, COMPASS) determined the axial charge a of the proton. In the simple 
valence-quark parton model, this can be identified with the quark spin and its observed 
suppression led to an intense experimental and theoretical search over two decades for the 
origin of the proton spin. In fact, as Veneziano was the first to understand [4], a does not 
measure spin in QCD itself and its suppression is related to OZI violations induced by the 
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U(1)a anomaly. 

In a series of papers, summarised in section 5, we have shown how a decouples from 
the real angular momentum sum rule for the proton (the form factors for this sum rule 
are given by generalised parton distributions (GPDs) which can be extracted from less 
inclusive measurements such as deeply-virtual Compton scattering) and is instead related 
to the gluon topological susceptibility [8, 9]. The experimentally observed suppression is a 
manifestation of topological charge screening in the QCD vacuum. In a 1994 paper with 
Narison [10], using QCD spectral sum rule methods, we were able to compute the slope 
of the topological susceptibility and give a quantitative prediction for a . Our prediction, 
a = 0.33, has within the last few months been spectacularly confirmed by the latest data 
on the deuteron structure function from the COMPASS and HERMES collaborations. 

Hopefully, this impressive new evidence for topological charge screening will provide 
fresh impetus to experimental 'spin' physics - first, to verify the real angular momentum 
sum rule by measuring the relevant GPDs, and second, to pursue the programme of target- 
fragmentation studies in semi-inclusive DIS at polarised ep colliders which we have proposed 
as a further test of our understanding of the g{ sum rule [11]. 

This review has been prepared in celebration of the 65th birthday of Gabriele Veneziano. 
I first met Gabriele when I came to Geneva as a CERN fellow in 1981. In fact, our first 
interaction was across a tennis court, in a regular Friday doubles match with Daniele Amati 
and Toine Van Proeyen. I like to think that in those days I could show Gabriele a thing 
or two about tennis - physics, of course, was a different matter. It has been my privilege 
through these ensuing 25 years to collaborate with one of the most brilliant and innova- 
tive physicists of our generation. But it has also been fun. As all his collaborators will 
testify, his good humour, generosity to younger colleagues, and enthusiasm in thinking out 
solutions to the deepest and most fundamental problems in particle physics and cosmology 
make working with Gabriele not only intellectually rewarding but hugely enjoyable. 

In his contribution to the 'Okubofest' in 1990 [12], Gabriele concluded an account of 
the relevance of the OZI rule to g\ by hoping that he had 'made Professor Okubo happy'. 
In turn, I hope that this review will make Gabriele happy: happy to recall how his original 
ideas on the U{\)a problem have grown into a quantitative description of anomalous QCD 
phenomenology, and happy at the prospect of new discoveries from a rich programme 
of experimental physics at future polarised colliders. It is my pleasure to join all the 
contributors to this volume in wishing him a happy birthday. 

2. The U(1)a anomaly and the topological susceptibility 

We begin by reviewing some essential features of the U(1)a anomaly, chiral Ward identities 
and the renormalisation group, placing particular emphasis on the role of the gluon topo- 
logical susceptibility. As we shall see, the anomaly provides the vital link between quark 
dynamics and gluon topology which is essential in understanding a range of phenomena in 
polarised QCD phenomenology. 
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2.1 Anomalous chiral Ward identities 

An anomaly arises when a symmetry which is present in the classical limit cannot be 
consistently imposed in a quantum field theory. The original example of an anomaly, and 
one which continues to have far-reaching implications for the phenomenology of QCD, is 
the famous Adler-Bell-Jackiw axial anomaly [13, 14, 15], which was first understood in 
its present form in 1969. In fact, calculations exhibiting what we now recognise as the 
anomaly had already been performed much earlier by Steinberger in his analysis of meson 
decays [16] and by Schwinger [17]. 

Anomalies manifest themselves in a number of ways. The original derivations of the 
axial anomaly involved the impossibility of simultaneously imposing conservation of both 
vector and axial currents due to regularisation issues in the AW triangle diagram in QED. 
More generally, they arise as anomalous contributions to the commutation relations in 
current algebra. A modern viewpoint, due to Fujikawa [18], sees anomalies as due to 
the non-invariance of the fermionic measure in the path integral under transformations 
corresponding to a symmetry of the classical Lagrangian. In this approach, the result of a 
chiral transformation q — > e iaaTa " ys q on the quark fields in the QCD generating functional 
W[V« 5 , V«, 9, S$, S a ] defined as 1 
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The terms in the square bracket are simply those arising from Noether's theorem, including 
soft breaking by the quark masses, with the addition of the anomaly involving the gluon 
topological charge density Q. Re-expressing the chiral variation of the elementary fields in 
terms of a variation with respect to the sources V^ 5 , V£, 9, Sf , S a then gives the functional 
form of the anomalous chiral Ward identities: 



d » Wv Zs -\fZn~f5aoW e - d abc m b W S i 



+fabcV b W V c + fabcV^Wyc + d abc S b W S c - d abc S b W S c = 



(2.3) 



x Our notation follows that of ref.[3]. The currents and pseudoscalar fields J° 5 , Q, <j)% together with the 
scalar <f> a are defined by 



J%5 = qi»lr,T a q 
4>t = qi 5 T a q 



J a — mo. 

M = TI^T q 

4> a = qT a q 



Q — —trGnvG^ 1 " 

07T 



where G M „ is the field strength for the gluon field. Here, T l = ^A 1 are flavour SU(nf) generators, and 
we include the singlet U(1)a generator T° = l/^/2n/ and let the index a = 0,i. With this normalisation, 
tvT a T b = i<5 a6 for all the generators T a . This accounts for the rather unconventional factor \J'lnf in 
the anomaly equation but has the advantage of giving a consistent normalisation to the full set of decay 
constants including the flavour singlets f° v and f° v . 

We will only need to consider fields where i corresponds to a generator in the Cartan sub-algebra, so that 
a = 3, 8, for rij = 3 quark flavours. We define d-symbols by {T a ,T b } = d a b c T c . For nj = 3, the explicit 
values are dooo = do33 = doss = ^330 = dsso = \/2/3, (fes = ^383 = —^888 = ^1/3. 
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where we have abbreviated functional derivatives as suffices. This is the key to all the 
results derived in this section. It makes precise the familiar statement of the anomaly as 



&>J% - y / 2n~fQ5 a0 - d abc m b ^ 5 



~ 



(2.4) 



The chiral Ward identities for two and higher-point Green functions are found by 
taking functional derivatives of eq.(2.3) with respect to the sources. The complete set of 
identities for two-point functions is given in our review [19]. As an example, we find 2 



3"<0|T J£ 5 <^|0> - y/^TfS^T Q 4\0) - d adc m d (0\T <f> c 5 4\0) - d abc {cf> c ) = (2.6) 



The anomaly breaks the original U(rif)L x U{nf)n chiral symmetry to SU(n,f)L x 
SU(nf)n x U(l)v/Z^ f and the quark condensate spontaneously breaks this further to the 
coset SU(nf)L>< SU(nf)/SU(rif)v- Goldstone's theorem follows immediately. In the chiral 
limit, there are (nj — 1) massless Nambu-Goldstone bosons, which acquire masses of order 
yfw, for non-zero quark mass. There is no flavour singlet Nambu-Goldstone boson since 
the corresponding current is anomalous. 

The zero-momentum Ward identities are especially important here, since they control 
the low-energy dynamics. With the assumption that there are no exactly massless particles 
coupling to the currents, we find 



Another key element of our analysis will be the chiral Ward identities for the effective 
action T[V^ 5 , V®, Q, <p a ], defined as the generating functional for vertices which are 1PI 
with respect to the set of fields Q, 4% and 4> a but not the currents J^ 5 , J^. This is achieved 
using the partial Legendre transform (or Zumino transform): 



d » W V" 6 S b 6 ~ 



y/2nfSaaW 0S b - M ac W S c S b - <S> ab = 



(2.5) 



which in more familiar notation reads 



y/^nJSaaWee + M ac W Si e = 



\/2nf8 a0 W 9S b + M ac W S c S b + $ ab = 



(2.7) 




2 We use the following 577(3) notation for the quark masses and condensates: 




and 



/(uu) \ 

{dd} = 2 ^ (<p a )T a 

\ (SS) ) 0.3,8 



where (<f> a ) is the VEV of <f> a = q T a q. It is also convenient to use the compact notation 



M ab = d acb m 



$ab = d abc {(f> C ) 
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The chiral Ward identities for T are 



d » Vv fe ~\/2nf S aoQ ~ d ahc m b (j)l 



+fabcVfr V c 5 + fabcV^Tyc - cW^I> + d abc <p C T^ = (2.9) 



Again, the zero-momentum identities for the two-point vertices play an important role: 

®acT<t>lQ - ^f2n~f5 a o = 

® ac T m -M ab = Q (2.10) 

These will be used in section 3 to construct an effective action which captures the low- 
energy dynamics of QCD in the pseudoscalar sector. 

2.2 Topological susceptibility 

The connection with topology arises through the identification of the gluon operator Q in 
the anomaly with a topological charge density. Q is a total divergence: 

Q = ^ tr G^CT = (2.11) 

where is the Chern-Simons current, 

K, = ^e, up MA u G^ - \gA»[AP, A"]) (2.12) 

Nevertheless, the integral over (Euclidean) spacetime of Q need not vanish. In fact, for 
gauge field configurations such as instantons which become pure gauge at infinity, 



(fx Q = neZ (2.13) 

where the integer n is the topological winding number, an element of the homotopy group 
ir 3 (SU(N c )). 

The form of the anomaly is then understood as follows. Under a chiral transformation, 
the fermion measure in the path integral transforms as (for one flavour) 

VqVq -► e -2i<*Jdx<phwi VqrDq = exp -2ia(n+-n_) VqVq (2.14) 

where (fi is a basis of eigenfunctions of the Dirac operator p in the background gauge 
field. The non-zero eigenvalues are chirality paired, so the Jacobian only depends on the 
difference (n+ — n_) of the positive and negative chirality zero modes of ]J). Finally, the 
index theorem relates the anomaly to the topological charge density: 

indp = n+-n- = J d A x Q (2.15) 
The topological susceptibility x{p 2 ) is defined as the two-point Green function of Q, 



viz. 



X {p 2 ) = ijdx e^(0|T Q(x) Q(0)|0) (2.16) 
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We are primarily concerned with the zero-momentum limit x(0) = Wgg(0). Combining 
eqs.(2.7) gives the crucial Ward identity satisfied by x(0): 



2n /X (0) = M 0a W S a Sl M m + (M<D) 00 (2.17) 

that is, 

J dx (0\T Q(x) Q(0)|0) = J dx m a m b (0\T $?(x) </> b 5 (0)\0) + m a (<p a ) (2.18) 



Determining exactly how this is satisfied in QCD is at the heart of the Witten-Veneziano 
approach to the U(1)a problem [20, 1]. 

The zero-momentum Ward identities allow us to write a precise form for the topological 
susceptibility in QCD in terms of just one unknown dynamical constant [21]. To derive 
this, recall that the matrix of two-point vertices is simply the inverse of the two-point Green 
function matrix, so in the pseudoscalar sector we have the following inversion formula: 

T QQ = -(Wee - Wes^Ws^rJWs^y 1 (2.19) 
Using the identities (2.7) and (2.17), this implies that at zero momentum 

F QQ = " ^fXiM^y 1 (2.20) 

and inverting this relation gives 

x = -r^ 1 " 2 ™/ r w( M$ )oo 1 )~ 1 (2-21) 

Finally, substituting for (M&)qq using the definitions above, we find the following impor- 
tant identity which determines the quark mass dependence of the topological susceptibility 
in QCD: 

*°>=-4- A ?4i>r <2 - 22) 

where we identify the non-perturbative coefficient A as Tqq. 

Notice immediately how this expression exposes the well-known result that x(0) van- 
ishes if any quark mass is set to zero. In section 3, we will see how it also clarifies the role 
of the l/N c expansion in the U(1)a problem. 

2.3 Renormalisation group 

The conserved current corresponding to a non-anomalous symmetry is not renormalised 
and has vanishing anomalous dimension. However, an anomalous current such as the 
flavour singlet axial current J^ 5 is renormalised. The composite operator renormalisation 
and mixing in the J° 5 , Q sector is as follows [22] : 

Qr = Qb~ -7=(1 - Z)d»J% B (2.23) 
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Notice the form of the mixing of the operator Q with d^J^ 5 under renormalisation. This 
ensures that the combination (<9^J° 5 — y / 2njQ) occurring in the U(1)a anomaly equa- 
tion is RG invariant. The chiral Ward identities therefore take precisely the same form 
expressed in terms of the bare or renormalised operators, making precise the notion of 
'non-renormalisation of the anomaly'. We may therefore interpret the above Ward identi- 
ties, which were derived in terms of the bare operators, as identities for the renormalised 
composite operators (and omit the suffix r for notational simplicity). 

The renormalisation group equation (RGE) for the generating functional W[V^ 5 , V®,9, S£, S a ] 
follows immediately from the definitions (2.23) of the renormalised composite operators. 
Including also a standard multiplicative renormalisation = Z m 1 for the pseudoscalar 
and scalar operators 4>§ and 4> a and denoting the anomalous dimensions corresponding to 
Z and Z^ by 7 and 7^ respectively, we find 3 

VW = 7(^5 - -j^d^Wvfr + 7</>(s 5 a Wsa + S a W s «) + • • • (2.24) 
where V = + 0% - 7m £ ? m q ^ 



V,8,S 5 ,S 



The RGEs for Green functions are found by functional differentiation of eq.(2.24) and 
can be simplified using the Ward identities. For example, for Wgg we find 

VW ee = 2>yWo e + 2 7 — ^=M ob W esi + . . . (2.25) 

V 2n f 



At zero momentum, we can then use the first identity in eq.(2.7) to prove that the topo- 
logical susceptibility %(0) is RG invariant, 

2? X (0) = (2.26) 

which is consistent with its explicit expression (2.22). 

A similar RGE holds for the effective action T[V^ 5 ,V^,Q,4>f,4>^\, which allows the 
scaling behaviour of the proper vertices involving Q and ^5 to be determined [23, 9, 24]. 
This reads 

vr = 7 (v;° 5 - -^=r Q ^)r y o 5 - 7^% + <p a r r ) + ... (2.27) 

An immediate consequence is that T>Tqq = at zero momentum, which ensures the com- 
patibility of (2.22) with the RG invariance of x(0)- 

2.4 l/N c , the topological expansion and OZI 

The final theoretical input into our analysis of the U(1)a problem and phenomenological 
implications of the anomaly concerns the range of dynamical approximation schemes asso- 
ciated with the large- N c limit. At various points we will refer either to the original large- N c 



3 The notation + . . . refers to additional terms which are required to produce the contact term contri- 
butions to the RGEs for n-point Green functions and vertices of composite operators. These are discussed 
fully in refs.[23, 9, 24], but will be omitted here for simplicity. They vanish at zero-momentum. 
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Figure 1: A typical Feynman diagram allowed in the large- N c limit. The dots on the quark loop 
represent external sources. 




Figure 2: Feynman diagrams allowed in the quenched approximation (left) or leading order in the 
topological expansion (right). 

expansion of 't Hooft [25], the topological expansion introduced by Veneziano [26] and the 
OZI limit [27, 28, 29]. A very clear summary of the distinction between them is given in 
Veneziano's 'Okubofest' review [12], which we follow here. 

In terms of Feynman diagrams, the leading order in the large N c , fixed rif ('t Hooft) 
limit is the most restrictive of these approximations, including only planar diagrams with 
sources on a single quark line and no further quark loops (Fig. 1). 

A better approximation to QCD is the quenched approximation familiar in lattice 
gauge theory. This is a small ny expansion at fixed iV c , i.e. excluding quark loops but 
allowing non-planar diagrams (Fig. 2). 

An alternative is the topological expansion, which allows any number of internal quark 
loops, but restricts to planar diagrams at leading order. Provided the sources remain 
attached to the same quark line, this corresponds to taking large N c at fixed nf/N c . This 
means that quarks and gluons are treated democratically and the order of approximation 
is determined solely by the topology of the diagrams (Fig. 2). 

Finally, the OZI approximation is a still closer match to full QCD with dynamical 
quarks than either the leading order quenched or topological expansions. Non-planar di- 
agrams and quark loops are retained, but diagrams in which the external sources are 
connected to different quark loops are still excluded (Fig. 3). This means that amplitudes 
which involve purely gluonic intermediate states are suppressed. This is the field-theoretic 
basis for the original empirical OZI rule. 

In each of these large- N c expansions, except the topological expansion where nf/N c 
is fixed, the U(1)a anomaly does not contribute at leading order. More precisely, the 
anomalous contribution (0\T Q <fi§\0} in the chiral Ward identity (2.6) is suppressed by 
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Figure 3: Feynman diagrams allowed (left) and forbidden (right) by the OZI rule. 



0(l/N c ) relative to the current term (0|T J° 5 0| 1 0) . This means that the flavour singlet 
current is conserved, Goldstone's theorem applies, and conventional PCAC methods can be 
used to understand the dynamics of the Green functions with a full set of (nj — 1) massless 
bosons in the chiral limit. Taking this as a starting point, we can then learn about the 
spectral decomposition of the actual QCD Green functions as we relax from the leading- 
order limits. In particular, this leads us to the famous Witten-Veneziano mass formula for 
the rf meson [20, 1]. 

The behaviour of the topological susceptibility at large N c is central to this analysis. 
It is clear from looking at planar diagrams that at leading order in 1/N C , x(0) in QCD 
coincides with the topological susceptibility x(0)|ym in the corresponding pure Yang-Mills 
theory. Referring now to the explicit expression (2.22) for x(0)> large- N c counting rules 
give A = 0(1) while (qq) = 0(N C ). It follows that for non-zero quark masses, 

X(0) = -A + 0{n f /N c ) (2.28) 

where A = Tqq is identified as — x(0)ym + 0(1 /N c ). On the other hand, if we consider the 
limit of x(0) for m q — > at finite N c , then we have 

X(0)U^ = (2.29) 

The 't Hooft large-iV c limit is therefore not smooth in QCD; the N c — > oo and m q — ► 
limits do not commute [20, 1, 21]. This is remedied in the topological expansion, where 
quark loops are retained and the 0(nf/N c ) contribution in eq.(2.28) allows the smooth 
chiral limit x(0) — ► even for large N c . 

3. 1 U(1)a without instantons' 

The U(1)a problem has a long history, pre-dating QCD itself, and has been an important 
stimulus to new theoretical ideas involving anomalies and gluon topology. 

At its simplest, the original 1 U(1)a problem' in current algebra is relatively straight- 
forwardly resolved by the existence of the anomalous contributions to the chiral Ward 
identities (anomalous commutators in current algebra) and the consequent absence of a 
ninth light Nambu-Goldstone boson in nj = 3 QCD. 4 However, a full resolution requires a 

4 The existence of a light flavour- singlet Nambu-Goldstone boson would produce a rapid off-shell varia- 
tion in the rj — > 3tt decay amplitude, in contradiction with the experimental data [30]. 
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much more detailed understanding of the dynamics of the pseudoscalar sector and the role 
of topological fluctuations in the anomalous Green functions. 

In this section, we review the analysis of the U(1)a problem presented by Veneziano in 
his seminal 1974 paper, 1 U(1)a without instantons' [l]. 5 As well as deriving the eponymous 
mass formula relating the rf mass to the topological susceptibility, the essential problem 
resolved in ref. [1] is how to describe the dynamics of the Green functions of the pseudoscalar 
operators in QCD in terms of a spectral decomposition compatible with the n/, N c , 6 and 
quark mass dependence imposed by the anomalous Ward identities. 

First, recall that in the absence of the anomaly, there will be light pseudoscalar mesons 
r/ a coupling derivatively to the currents with decay constants f aa , i.e. (0| J^\i] a ) = ip^f aa - 
(We use the notation rf 1 to denote the physical mesons vr°,ry and rf, while the SU(3) 
index a = 3,8,0.) The mass matrix \j? a ^ satisfies the Dashen, Gell-Mann-Oakes-Renner 
(DGMOR) relation [35, 36] 

r»lpf Pb = ~ (M$U (no anomaly) (3.1) 

The consequences of the anomaly are determined by the interaction of the pseudoscalar 
fields 05 with the topological charge density Q and the subsequent mixing. This gives 
rise to an additional contribution to the masses. Moreover, we can no longer identify 
the flavour singlet decay constant by the coupling to J^ 5 since this is not RG invariant. 
Instead, the physical decay constants f aa are defined in terms of the couplings of the rf to 
the pseudoscalar fields through the relation f aa {Q\4> b 5 \r] a ) = d a b c (<fi c )- This coincides with 
the usual definition except in the flavour singlet case. 

The most transparent way to describe how all this works is to use an effective action 
r[Q,(/)|] constructed to satisfy the anomalous chiral Ward identities. It is important to 
emphasise from the outset that this is an effective action in the sense of section 2.1, i.e. the 
generating functional for vertices which are 1PI with respect to the set of fields Q,4>§ 
only. The choice of fields is designed to capture the degrees of freedom essential for the 
dynamics. 6 A different choice (or linear combination) redefines the physical meaning of 
the vertices, so it is important that the final choice of fields in T results in vertices which 
are most directly related to physical couplings. 

The simplest effective action consistent with the anomalous Ward identities and the 
renormalisation group is 

T[Q,r 5 ] = Jdx (JjQ 2 + Qi^Soa-Bad 2 )^ 1 ^ 

+ ^acHiM^U-C^d 2 )^ 1 ^ (3.2) 
5 For reviews of the instanton approach to the resolution of the U(1)a problem, see e.g. refs.[31, 32, 33, 

34]. 

6 Note especially the frequently misunderstood point that the choice of fields in F is not required to be in 
any sense a complete set, nor does the restriction to a given set of fields constitute an approximation. Before 
imposing dynamical simplifications, the identities derived from F are exact - increasing the set of basis fields 
simply changes the definitions of the 1PI vertices. The effective action considered here is therefore different 
from the non-linear chiral Lagrangians incorporating the large-iV c approach to the pseudoscalar mesons 
constructed by a number of groups. See, for example, refs. [37, 21, 38, 39, 40, 41, 42]. 
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The constants C ab and B a are related to T va V b and Yy a q respectively. The inclusion of 
the term with B a is unusual, but is required for consistency with the RGEs derived from 
(2.27) beyond zero momentum. 

This form of T[Q, </>g] encodes three key dynamical assumptions: 

• Pole dominance. We assume that the Green functions are dominated by the contribu- 
tion of single-particle poles associated with the pseudoscalar mesons including the flavour 
singlet. This extends the usual PCAC assumption to the singlet sector. 

• Smoothness. We assume that pole-free dynamical quantities such as the decay constants 
and couplings (1PI vertices) are only weakly momentum-dependent in the range from p = 
to their on-shell values. This allows us to impose relations derived from the zero-momentum 
Ward identities, provided this is compatible with the renormalisation group. 

• Topology. There must exist topologically non-trivial fluctuations which can give a non- 
vanishing value to x(0)|ym hi pure gluodynamics. This is required to give the non- vanishing 
coefficient in the all-important ^jQ 2 term in r[Q, $5]. Notice that we do not require 
a kinetic term for Q, which would be associated with a (presumed heavy) pseudoscalar 
glueball. 

The second derivatives of T[Q, 4>f] are 

Tqq r^*\ / A' 1 {^hT f 5 0d + B d p 2 )^ 

%q T m j UacHT^co + ^P 2 ) ^{{M®U + C cd p 2 )$>- d l 

(3-3) 

The corresponding Green functions (composite operator propagators) are given by inver- 
sion: 



Wee W 0S b \ _( T QQ F Q4> b x 



(3.4) 



and we find, to leading order in p 



,2 



Wee = -A A' 1 

W eS b = W s > e ~ y^AA^^b 

W S a Sl = -<5> ac Aj<S> db (3.5) 

where 

A = 1 — (2n f AS a0 S 0b + ^/2n~ f A(5 a0 B b + B a S 0b )p 2 ^j (M$ + Cp 2 )^ 1 (3.6) 

and 

A a6 = (C ab - y f2n~ f A(5 a0 B b + B a 5 0b ))p 2 + (M<S>) ab - 2n f A 5 a0 5 0b (3.7) 

In this form, however, the propagator matrix is not diagonal and the operators are not 
normalised so as to couple with unit decay constants to the physical states. It is therefore 
convenient to make a change of variables in T so that it is written in terms of operators 
which are more closely identified with the physical states. We do this is in two stages, since 
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the intermediate stage allows us to make direct contact with the discussion in ref.[l] and 
will play an important role in some of the phenomenological applications considered later. 

First, we define rescaled fields fj a whose kinetic terms, before mixing with Q, are 
canonically normalised. That is, we set 

with the 'decay constants' f aa defined such that ^r^p \ p =q = 5 a p. This implies 

d 

: Woa ob 

p=0 

where D a = y / 2rif5 a oQ + M ab (f>\ is the divergence of the current J° 5 . In the chiral limit, 
this reduces in the flavour singlet sector to 

(/7 T )oo = 4^X{P 2 ) „ = X'(0) (3.10) 



(ff 1 U = C ab = jpWgafr n (3.9) 



dp 



p=0 



a result which plays a vital role in understanding the 'proton spin' problem. Notice however 
that the f aa are not RG invariant: in fact, Vf aa = r yS a of aa . The effective action r[Q,rj a ] 
is: 

r[Q,?f] = Jdx (J^q 2 + Q{^M f 5» a -B a d 2 )u- i rr 

+ \fl a (-d 2 + f- 1T M<S>r%^) (3.11) 

In this form, the fj a are the canonically normalised fields corresponding to the 'would-be 
Nambu-Goldstone bosons' in the absence of the anomaly, before they acquire an additional 
anomaly-induced mass. In the framework of the large-iV c or OZI approximations, they 
would correspond to true Nambu-Goldstone bosons. The singlet fp is what we have there- 
fore referred to in our previous papers as the 'OZI boson' t{q ZI . As we see later, the naive 
current algebra relations hold when expressed in terms of the fj a and f aa , though these do 
not correspond to physical states or decay constants. 

The physical particle masses are identified with the poles in the two-point Green func- 
tions (3.5). We immediately see that due to mixing with the topological charge density Q, 
the physical pseudoscalar meson mass m 2 al3 is shifted from its original value. From the pole 
in eq.(3.7), we immediately find 

f aa rnlpfP b = -{M^>) ab + 271/^00*60 (3-12) 

where we identify the physical, RG-invariant decay constants as 

(ff T )ab = (ff T )ab-V2^~fA(5 a0 B b + B a 5 0b ) (3.13) 

Eq.(3.12) is the key result. It generalises the original DGMOR relations (3.1) to the 
flavour-singlet sector with the anomaly properly incorporated and the renormalisation 
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group constraints satisfied. It represents a generalisation of the Witten-Veneziano mass 
formula which makes no direct reference to large-iV c arguments but depends only on the 
three dynamical assumptions stated above [2]. 

With this clarification of the distinction between the physical decay constants f aa 
and the RG non-invariant f aa , we can rewrite eq.(3.6) for the topological susceptibility 
x {p 2 ) = W ee (p 2 ) as 

X (p 2 ) = -^[l-tr((/7 T -//V + 2n / Al 00 )(//V + M$)- 1 ] _1 (3.14) 

It is clear that in the zero-momentum limit, this expression successfully reproduces eq.(2.22) 
for x(0). For one flavour, the formula simplifies to 

X(p 2 ) = -A(ff T p 2 + M$)[ff T p 2 + M$ + 2n f AY 1 (n f = 1) (3.15) 

showing clearly the pole at the shifted mass m 2 of eq.(3.12). The occurrence of both f aa 
and f aa in these expressions allows them to satisfy the RGE (2.25) for the topological 
susceptibility, which requires Vx{p 2 ) = 0(p 2 ). 

The second stage is to make a change of variable which diagonalises the propagator 
matrix, so as to give the most direct possible relation between the operators and the 
physical states. Choosing 



G = Q-WeszW-}^ ~ Q + ^M f A^4>\ 

5 5 

V a = f Taa Kb4 (3-16) 
defines the effective action T[G,rj a ] as 

r[G,r, a ] = fdx {^G 2 + lrf(-d 2 -m 2 ) Q/3 /) (3.17) 

with m 2 a p given by eq.(3.12). The corresponding propagators are 

(0|T G G\0) = -A 

(0\T r] a rflO) = 9 ~ l 9 5 a(B (3.18) 
p z - m^ a 

where with no loss of generality we have taken m 2 ^ diagonal. 

Notice also that the states mix in the complementary way to the operators. In partic- 
ular, the mixing for the states corresponding to eq.(3.16) for the fields G and rf is 

\G) = \Q) 

\v a ) = (.T'n^a^l) - v^fA8 a0 \Q)) (3.19) 

In this sense, we see that we can regard the physical rf (and, with SU(3) breaking, the rj) 
as an admixture of quark and gluon components, while the unphysical state \G) is purely 
gluonic. 



- 14 - 



An immediate corollary is the following relation, which we will use repeatedly in de- 
riving alternative forms of the current algebra identities for the pseudoscalar mesons: 

The formulation in terms of T[G, r] a ] is exactly what we need to construct a simple 
1 U(1)a PCAC with which to interpret the low-energy phenomenology of the pseudoscalar 
mesons. We turn to this in the next section. 

Here, we focus on the intermediate formulation T[Q, fj a ] in order to describe Veneziano's 
analysis of the U (I) a problem in the framework of the large- N c and topological expansions. 
The starting point is the anomalous Ward identity (2.18) for the topological susceptibility: 



n) 



dx (0\T Q(x) Q(0)\Q) = j dx m a m b (0\T <j>%{x) </>|(0)|0) + m a ((f> a ) (3.21) 



The essential problem is how to understand this relation in terms of a spectral decompo- 
sition in the context of the 1/N C expansion. 

Assuming that x(0)ym = —A + 0(1/N C ) is non-vanishing at 0(1), the l.h.s. is 0(rij) 
in leading order in 1/N C . On the other hand, the r.h.s. includes the condensate term 
of 0(rifN c m). To resolve this apparent paradox, we have to go beyond leading order 
in 1/N C and consider the quark loop contributions which are included in the topological 
expansion. Although these are formally suppressed by powers of (rif/N c ), they contain 
light intermediate states which can enhance the order of the Green function. As we have 
seen above, these light states are just the 'OZI bosons' \i) a ) with masses //^ of 0(n/m). 
Inserting these intermediate states, we therefore find that: 

X(P 2 ) = x(p 2 )\YM-(0\Qm ( 2 \ ) (f\Q\0) + ... (3.22) 

\P A* )ct/3 

where the coupling (0\Q\rj a ) is 0{^nj /N c ). 

Approximating x(p 2 )ym ~ —A (a low-momentum smoothness assumption) and (0\Q\r) a ) 
y/2nfA(f~ 1 ) 0a , then summing the series of intermediate state contributions, we find 

X{P 2 ) ^ — r (3-23) 

l-2n / A(/(p2- M 2)/r) oo 

This expression reproduces eq.(13) of ref.fl]. Clearly, it is dominated by the physical 
pseudoscalar pole with anomaly- induced mass given by eq.(3.12). It does not completely 
recover our more precise expression (3.14) because of the approximation for the coupling of 
Q to the \f) a ), which misses the subtleties related to the introduction of B a in the effective 
action T[Q,fj a ] and the distinction of f aa and f aa . These are effects of higher order in 
l/N c but, as we have seen, they are necessary to establish full RG consistency and will 
prove to be important for phenomenology. 

To see how a term with the 0{rifN c m) dependence of the condensate can arise in 
n /X(0), notice from eq.(3.12) that the physical pseudoscalar mass squared m^ a has two 
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contributions, the first of 0(m) from the conventional quark mass term and the new, 
anomaly- induced contribution of 0(nf/N c ). If we are in a regime where the anomaly 
contribution dominates (m < A.qct>/N c ), then it follows that the above expression for x(0) 
indeed becomes of 0(n^ l N c m). 

The original Witten-Veneziano mass formula for the rj is the large- iV c limit of eq.(3.12). 
In the chiral limit there is no flavour mixing and the singlet mass is given by 



rn 



( -2n f A = -jJ- x (0) YM + O((n f /N c y) (3.24) 

This formula provided the first link between the rj mass and gluon topology. For an 
alternative recent derivation in the context of a nf/N c expansion, see also ref.[43]. 

What we learn from all this is that the Green functions in the anomalous chiral Ward 
identities admit a consistent spectral decomposition in terms of a full set of (rij — 1) 
pseudoscalar mesons, provided they satisfy the generalised DGMOR mass formula (3.12) 
including the all-important anomaly term. The presence of these light poles can enhance 
the apparent order of the Green functions, as is familiar with Nambu-Goldstone bosons, 
and the anomaly- induced 0(rif/N c ) contribution to m 2 a is critical in ensuring complete 
consistency with the Ward identities. 

Similar considerations apply to the resolution of apparent paradoxes in the ^-dependence 
of some Green functions. For example [1], we can show from the anomalous Ward identities 
that the condensate satisfies 

Y,™M)\e = m a {^ a ) = co S (6/n f ) m a (^)| e=0 (3.25) 



This implies 



d 2 



■m a (r)\e=o = -m a Jdxjdy (0\T Q(x) Q(y) a (O)|O) 



= -\ m ^^)\ e ^ (3.26) 

Here, the Green function is superficially of 0(n,f/N c ) while the r.h.s. is 0(N c /nf). The 
resolution is simply that it contains pseudoscalar intermediate states contributing two light 
poles with m 2 ~ 0(rif/N c ). So once again we see how the spectral decomposition in terms 
of the full set of pseudoscalar mesons, including the flavour singlet, ensures consistency 
with the anomalous Ward identities. 



4. Pseudoscalar mesons 

This theoretical analysis provides the basis for an extension of the conventional PCAC 
or chiral Lagrangian description of the phenomenology of the pseudoscalar mesons to the 
flavour singlet sector. In this section 7 we describe the role of the U(1)a anomaly in the 
radiative decays of 7r°,r/ and rj and derive the U(1)a Goldberger-Treiman relation, first 
proposed by Veneziano as a resolution of the 'proton spin' problem. 

7 This section is based on the presentation in ref. [3], where we extend and update our original work 
[2, 44] to include a detailed comparison with experimental data. 
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4.1 U(1)a Dashen, Gell-Mann Oakes Renner relations 

The extension of the DGMOR relations to the U(1)a sector follows from the application 
of the three key dynamical assumptions used above (viz. pole dominance by the nonet of 
pseudoscalar mesons, smoothness of decay constants and couplings over the range from 
zero to on-shell momentum, and the existence of topologically non-trivial gluon dynamics) 
to the anomalous chiral Ward identities. 

The fundamental U(1) A DGMOR relation 

f aa m 2 al3 f T ? b = -M ac $> cb + 2n f A5 a0 5 b0 (4.1) 

has been derived above in the course of the general discussion of the U(1)a problem. In 
order to make this section self-contained, we give a brief and direct derivation here. 

Recall that the physical meson fields are given as r/ a = f Taa ^~ b l (f)\, with the decay 
constants defined so that the propagator W SaS & = — l/(p 2 — m^) a p. It follows immediately 
that at zero momentum, 

f aa <pf Pb = *ac(W S5 S 5 & db (4.2) 

Using the chiral Ward identities of section 2 together with the identification (2.21) of the 
topological susceptibility, we can then show 

QaciWstsJzt**, = ($M) ac (MW S5 s 5 M)J(M<!>U 

= {M$)J-{M$) + 2n /X (0)l 00 r 1 (M$U 

V / cd 

= -(M*U + 2n f T Q 1 Q S a0 S b0 (4.3) 

proving the result (4.1). 

Expanding this out, and assuming the mixed decay constants / 07r , f 8lT , f 3 ' n , f 3 ' n ' are all 
negligible, we have 

(/ 0r 2 m2, + (/°") 2 m2 = -l(m u {uu)+m d {dd)+m s {ss)) + QA (4.4) 

fOv' fn' m >2 + f°vfv m 2 = _ v ( mu (^ n ) + md (dd) - 2m s (ss}) (4.5) 

(/V)2 m « + (/*?)2 m 2 = -±( mu (uu)+m d (dd)+4m s (ss)) (4.6) 

flml = - {m u (uu} + m d {dd)) (4.7) 

and we can add the standard DGMOR relation for the K + , 

fK m K = - (m u (uu) + m s (ss)) (4.8) 

We emphasise that these formulae, as well as the radiative decay and U(1) A Goldberger- 
Treiman relations derived below, do not depend at all on the 1/N C expansion. In particular, 
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the constant A appearing in the flavour singlet formula is defined as the non-perturbative 
parameter determining the topological susceptibility x(0) in QCD according to the exact 
identity (2.22). As explained above, large- N c ideas do indeed provide a rationale for ex- 
tending the familiar PCAC assumptions of pole dominance and smoothness to the flavour 
singlet channel, but these assumptions can be tested independently against experimental 
data. 

The most useful form of these relations for phenomenology is to assume exact SU (2) 
flavour symmetry and eliminate the quark masses and condensates in favour of f ni fx, m 2 
and m 2 K in the DGMOR relations for the r\ and rj '. This gives 

(/ (tf)2 m 2 +(/ o, )2m 2 = l(fiml + 2fj c m 2 K ) + 6A (4.9) 

WH 2 + /°v 8 X = (4.io) 

(/ 8V)2 m ,2 + (/ 8, ) 2 m 2 = - l - { flml-AfW K ) (4.11) 



We can also now clarify the precise relation of these results to the Witten-Veneziano 
formula for the mass of the rj' in its non-vanishing quark mass form, viz. 

m 2 ,+m 2 -2m 2 K = - — X (0)\ YM (4.12) 

Of course, only the m 2 , term on the l.h.s. is present in the chiral limit. Substituting in 
the explicit values for the masses in this formula gives a prediction [1] for the topological 
susceptibility, x(0) ~ —(180 MeV) 4 , which as we see below is remarkably close to the 
subsequently calculated lattice result. 

If we now add the DGMOR formulae (4.9) and (4.11), we find 

if*?™ 2 , + (n 2 m 2 + (D 2 ^ + (f^'frn 2 , - 2f 2 K m 2 K = 6A (4.13) 

which we repeat is valid to all orders in 1/N C . To reduce this to its Witten-Veneziano 
approximation, we impose the large-iV c limit to approximate the QCD topological charge 
parameter A with — x(0)|ym as explained in section 2.4. We then set the 'mixed' decay 
constants f° v and f 8v ' to zero and all the remaining decay constants f ' n ' , f 8v and fx 
equal to f w . With these approximations, we recover eq.(4.12). Eventually, after we have 
found explicit experimental values for all these quantities, we will be able to demonstrate 
quantitatively how good an approximation the large--/V c Witten-Veneziano formula is to 
the generalised U(l) A DGMOR relation in full QCD. 

4.2 Radiative decay formulae for vr°, r/, 77' — > 77 

Radiative decays of the pseudoscalar mesons are of particular interest as they are controlled 
by the electromagnetic ^7(1)^ anomaly, 

, rv ~ 

d"J^-M ab( /> b 5 -^2ir f Q6 a0 -a a em —F^F^ = (4.14) 
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where F^ v = d^Ay — dyA^ is the usual electromagnetic field strength and the anomaly coeffi- 
cients cig m are determined by the quark charges. The generating functional r[V^ 5 , V£, Q, <j>^, 4> a , A^] 
of 1PI vertices including the photon satisfies the Ward identity 



^ry« -^T f 5 a0 Q - a a em —F^F^ - d abc m b 4>% 

+f abc Vfr V c s + f abc V* 5 T V c - d abc ^T^ b + dabctT^ = (4.15) 

To derive the radiative decay formulae, we first differentiate this identity with respect 
to the photon field A^. This gives 

(X 

iPvFv\A*AP + ^ab^^AXAf = ~ a cm~ Z»v\ P K k 2 ( 4 - 16 ) 

where ki,k 2 are the momenta of the two photons. Notice that the mass term does not 
contribute directly to this formula. Prom its definition as 1PI w.r.t. the pseudoscalar 
fields, the vertex T'v%a x ap does not have a pole at p 2 = 0, even in the massless limit, so 
we find simply 

^a,b^^b A \ A p\ p=Q = ~ a tm~ e ^Xpk 1 k2 (4-17) 

The radiative couplings g^j-y for the physical mesons i] a = ir°,r],r]' are defined as 
usual from the decay amplitude (ttI 7 ?")- With the PC AC assumptions already discussed, 
they can be identified with the 1PI vertices as follows: 

<77lO = " ig^ry e^ Xp k^e\k 1 )eP(k 2 ) = iV r}aAxAp e x {k 1 )e p {k 2 ) (4.18) 

Re-expressing eq.(4.17) in terms of the canonically normalised 'OZI bosons' r) a , we therefore 
have the first form of the decay formula, 

r^ 77 = cC^ (4.19) 

Then, rewriting this in terms of the physical pseudoscalar couplings g^^ and decay 
constants according to the relation (3.20) gives the final form for the generalised U(1) A 
PCAC formula describing radiative pseudoscalar decays, incorporating both the electro- 
magnetic and colour anomalies: 



/ a %«77 + y^fAg Gll S a0 = a° m - (4.20) 
Expanding this formula, we have 

/""Wrr + /°%77 + yfcA 9Grr = a° em - (4.21) 



7T 



/ ^^'77 + / ^drfyy — a cm ^ (4.22) 
3 ^ 

/7r5 f 7r 7 7 = "cm" (4.23) 



where a° em = ^N c , a 8 cm = ^N c and a 3 cm = ±N C . 
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The new element in the flavour singlet decay formula is the gluonic coupling parameter 
gc-yy- It takes account of the fact that because of the anomaly- induced mixing with the 
gluon topological density Q, the physical rf is not a true Nambu-Goldstone boson so the 
naive PCAC formulae must be modified. 5G77 is not a physical coupling and must be 
regarded as an extra parameter to be fitted to data, although in view of the identifications 
in eq.(3.19) it may reasonably be thought of as the coupling of the photons to the gluonic 
component of the rj . 

The renormalisation group properties of these relations are readily derived from the 
RGE (2.27) for T. In the 'OZI boson' form, the unphysical coupling g^a 77 satisfies the 
complementary RGE to the decay constant f aa so the combination is RG invariant: 

Vf aa = jS a0 f aa V(f aa g r ^) = (4.24) 

In contrast, all the decay constants and couplings in the relation (4.20) can be shown to 
be separately RG invariant, including the gluonic coupling gcyy [24, 44]. 

4.3 The renormalisation group, OZI and 1/N C : a conjecture 

Although these U(1)a PCAC relations have been derived purely on the basis of the pole 
dominance and smoothness assumptions, we will nevertheless find it useful in practical 
applications to exploit their OZI or large- N c behaviour, in conjunction with the renormal- 
isation group. 

The basic idea is that violations of the OZI rule, or equivalently anomalous large-iV c 
behaviour, are generally related to the existence of the U{1)a anomaly. Moreover, we 
can identify the quantities which will be particularly sensitive to the anomaly as those 
which have RGEs involving the anomalous dimension 7. We therefore conjecture that 
the dependence of Green functions and 1PI vertices on 7 will be an important guide in 
identifying propagators and couplings which are likely to show violations of the OZI rule 
and those for which the OZI (or large- N c ) limit should be a good approximation [9, 24]. 

As an example, the large-iV c order of the quantities in the flavour singlet decay rela- 
tion (4.21) is as follows: f aa = 0(\/W c ) for all the decay constants, ^« 77 = 0(y/Nc), 
<7G77 = 0(1), a" m = 0(N C ) and the topological susceptibility parameter A = 0(1). The 
renormalisation group behaviour is especially simple, with both the meson and gluonic 
couplings <7r7« 77 and gcy^ as well as the decay constants being RG invariant. Putting this 
together, we find that all the terms in the decay formula are of 0(N C ) except the anoma- 
lous contribution AgQ^ which is 0(1). Since it is RG invariant and independent of the 
anomalous dimension 7, we conjecture that it is a quantity for which the OZI (or large- N c ) 
approximation should be reliable so we expect it to be numerically small compared with 
the other contributions. In the next section, we test this against experiment. 

As we shall see later, this conjecture has far-reaching implications for a range of predic- 
tions related to the anomaly, particularly in the interpretation of the U(\)a Goldberger- 
Treiman relation and associated ideas on the first moment sum rules for g\ and gj in 
deep-inelastic scattering. 
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4.4 Phenomenology 

After all this theoretical development, we finally turn to experiment and use the data on the 
radiative decays rj, rj — > 77 to deduce values for the pseudoscalar meson decay constants 
f° v ', f° v , f 8r> ' and f 8r ) from the set of decay formulae (4.21), (4.22) and U(1) A DGMOR 
relations (4.9)-(4.11). We will also find the value of the unphysical coupling parameter 
<7G 77 and test the realisation of the l/N c expansion in real QCD. 
The two-photon decay widths are given by 

3 

rtffa) -» 77) = ^K'(„) 77 | 2 (4.25) 

The current experimental data, quoted in the Particle Data Group tables [45], are 

r(n' -► 77) = 4.28 ± 0.19 KeV (4.26) 

which is dominated by the 1998 L3 data [46] on the two-photon formation of the rf in 
e+e~ — > e + e~7r + 7r~7, and 

r(r/ -► 77) = 0.510 ± 0.026 KeV (4.27) 

which arises principally from the 1988 Crystal Ball [47] and 1990 ASP [48] results on 
e + e~ — > e + e~7]. From this data, we deduce the following results for the couplings gy 77 

grfvy = 0.031 ± 0.001 GeV" 1 (4.28) 

and 

grrr, = 0.025 ± 0.001 GeV" 1 (4.29) 

which may be compared with g nll = 0.024 ± 0.001 GeV. 
We also require the pseudoscalar meson masses: 

rry = 957.78 ± 0.14 MeV m v = 547.30 ± 0.12 MeV 

m K = 493.68 ± 0.02 MeV = 139.57 ± 0.00 MeV (4.30) 

and the decay constants f v and fx- These are defined in the standard way, so we take the 
following values (in our normalisations) from the PDG [45]: 

f K = 113.00 ±1.03 MeV U = 92.42 ± 0.26 MeV (4.31) 

giving f K /f n = 1.223 ± 0.012. 

The octet decay constants f 8ri and f 8v are obtained from eqs.(4.11) and (4.22). This 
leaves three remaining equations which determine the singlet decay constants f° v ' , f° v and 
the gluonic coupling gc~ n in terms of the QCD topological susceptibility parameter A. 
This dependence is plotted in Figs. 4 and 5. 

To make a definite prediction, we need a theoretical input value for the topological 
susceptibility. In time, lattice calculations in full QCD with dynamical fermions should 
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Figure 4: The decay constants f° n and f° v as functions of the non-perturbative parameter A = 
(x MeV) 4 which determines the topological susceptibility in QCD. 
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Figure 5: This shows the relative sizes of the contributions to the flavour singlet radiative decay 
formula (4.21) expressed as functions of the topological susceptibility parameter A — (x MeV) 4 . 
The dotted (black) line denotes ^^^f 21 - The dominant contribution comes from the term f° v g^yy, 
denoted by the long-dashed (green) line, while the short-dashed (blue) line denotes f^g^yy- The 
contribution from the gluonic coupling, ^6^3077, is shown by the solid (red) line. 

be able to determine the parameter A. For the moment, however, only the topological 
susceptibility in pure Yang- Mills theory is known accurately. The most recent value [49] is 

x (0)\ym = - (191 ±5 MeV) 4 = - (1.33 ± 0.14) x 1(T 3 GeV 4 (4.32) 

This supersedes the original value x(^)\ym — —(180 MeV) 4 obtained some time ago [50]. 
Similar estimates are also obtained using QCD spectral sum rule methods [51]. At this 
point, therefore, we have to make an approximation and so we assume that the 0(1/N C ) 
corrections in the identification 

A = x(0)ym + 0(1/N C ) (4.33) 

are numerically small. With this provisional input for A, we can then determine the full 
set of decay constants: 

/ 0r? ' = 104.2 ± 4.0 MeV f° v = 22.8 ± 5.7 MeV 

fn' = -36.1 ±1.2 MeV f v = 98.4 ± 1.4 MeV (4.34) 

and 

5G77 = - 0.001 ± 0.072 GeV" 4 (4.35) 
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It is striking how close both the diagonal decay constants f 0r1 ' and f 8v are to /„-. Pre- 
dictably, the off-diagonal ones / 0r? ' and f Sv ' are strongly suppressed. 

It is also useful to quote these results in the two-angle parametrisation normally used 
in phenomenology. Defining, 



we find 



f° v ' f° v \ //ocos#o -/osin^o 

f v ' f Sr 'J \f 8 sme 8 f 8 cos8 8 



fo = 106.6 ± 4.2 MeV / 8 = 104.8 ± 1.3 MeV 



(4.36) 



that is 



9 = - 12.3 ± 3.0 deg 6 8 = - 20.1 ± 0.7 deg (4.37) 
j- = 1.15 ±0.05 ^ = 1.13 ±0.02 (4.38) 

Jir Jit 

Given these results, we can now investigate how closely our expectations based on OZI 
or 1/N C reasoning are actually realised by the experimental data. With the input value 
(4.32) for A, the numerical magnitudes and 1/N C orders of the terms in the flavour singlet 
decay relation are as follows (see Fig. 5): 



f^g^yr [N c \ 3.23] + /°% 77 [N c ; 0.57] + V6Ag Gll [1; - 0.005 ± 0.23] 
= a° m ^ [N c ; 3.79] (4.39) 

7T 

The important point is that the gluonic contribution gG-yy, which is suppressed by a power 
of 1/N C compared to the others, is also experimentally small. The near- vanishing for 
the chosen value of A is presumably a coincidence, but we see from Fig. 5 that across 
a reasonable range of values of the topological susceptibility it is still contributing no 
more than around 10%, in line with our expectations for a RG-invariant, OZI-suppressed 
quantity. 

It is also interesting to see how the 1/N C approximation is realised in the U(1)a DG- 
MOR generalisation (4.13) of the Witten-Veneziano formula (4.12). Here we find 

(/°"') 2 ^ [iV c ; 9.96] + (f 0r >) 2 m 2 [N c ; 0.15] + (f T >') 2 m' 2 [N c ; 1.19] 
+ (fvfml [N c ; 2.90] -2f 2 K m 2 K [N c - -6.22] = 6A [1; 7.98] 

(4.40) 

This confirms the picture that the anomaly-induced contribution of 0(1/N C ) to m 2 /, which 
gives a sub-leading 0(1) effect in (f 0ll ') 2 m 2 , , is in fact numerically dominant and matched 
by the O(l) topological susceptibility term 6A Away from the chiral limit, the conven- 
tional non-anomalous terms are all of 0(N C ) and balance as expected. The surprising 
numerical accuracy of the Witten-Veneziano formula (2.18) is seen to be in part due to a 
cancellation between the underestimates of f 8r >' (taken to be 0) and fx (set equal to f n ). 
This emphasises, however, that great care must be taken in using the formal order in the 
l/N c expansion as a guide to the numerical importance of a physical quantity, especially 
in the U{V)a channel. 
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Nevertheless, the fact that the RG-invariant, OZI-suppressed coupling is exper- 
imentally small is a very encouraging result. It increases our confidence that we are able 
to identify quantities where the OZI, or leading l/N c , approximation is likely to be nu- 
merically good. It also shows that <7g 77 gives a contribution to the decay formula which is 
entirely consistent with its picturesque interpretation as the coupling of the photons to the 
anomaly- induced gluonic component of the rf . A posteriori, the fact that its contribution is 
at most 10% explains the general success of previous theoretically inconsistent phenomeno- 
logical parametrisations of 7/ decays in which the naive current algebra formulae omitting 
the gluonic term are used. 

However, while the flavour singlet decay formula is well-defined and theoretically con- 
sistent, it is necessarily non-predictive. To be genuinely useful, we would need to find 
another process in which the same coupling enters. The problem here is that, unlike the 
decay constants which are universal, the coupling is process-specific just like g,/ 77 
or gryy^. There are of course many other processes to which our methods may be applied 
such as r]'(r]) — > V7, where V is a flavour singlet vector meson p,u,4>, or r/'(r/) — > 7r + 7r~7. 
The required flavour singlet formulae may readily be written down, generalising the naive 
PCAC formulae. However, each will introduce its own gluonic coupling, such as gcv^- 
Although strict predictivity is lost, our experience with the two-photon decays suggests 
that these extra couplings will give relatively small, at most O(10 — 20%), contributions 
if like gcyy they can be identified as RG invariant and 1/N C suppressed. This observation 
restores at least a reasonable degree of predictivity to the use of PCAC methods in the 
U(1)a sector. 

4.5 U(1)a Goldberger-Treiman relation 

A further classic application of PCAC is to the pseudoscalar couplings of the nucleon. For 
the pion, the relation between the axial-vector form factor of the nucleon and the pion- 
nucleon coupling g n NN is the famous Goldberger-Treiman relation. Here, we present its 
generalisation to the flavour singlet sector, which involves the anomaly and gluon topol- 
ogy. This U{\)a Goldberger-Treiman relation was first proposed by Veneziano [4] in an 
investigation of the 'proton-spin' problem and further developed in refs.[8, 9, 52, 3]. 
The axial- vector form factors are defined from 

(N\J- 5 \N) = 2m N (G\(p 2 )s^ + G a P {p 2 )p.s Pll ) (4.41) 

where = u^^^uj^m^ is the covariant spin vector. In the absence of a massless pseu- 
doscalar, only the form factors G A (0) contribute at zero momentum. 

Expressing the matrix element in terms of the 1PI vertices derived from the generating 
functional r[V^ 5 , V®, Q, cf) a ], including spectator fields N, N for the nucleon, we have 

(N\J^\ N ) = u( T vrNN + W vre T Qm + W vrs f mN y (4.42) 

Note that this expansion relies on the specific definition(2.8) of T as a partial Legendre 
transform. 
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We also need the following relation, valid for all momenta, which is derived directly 
from the fundamental anomalous chiral Ward identity (2.9) for T: 

d V T V« 5 NN = -^ab^^NN (4.43) 

Now, taking the divergence of eq.(4.42), using this Ward identity and then 8 taking the 
zero- momentum limit, noting that the propagators vanish at zero momentum since there 
is no massless pseudoscalar, gives 

2m N G a A {0) uj 5 u = iu <S>abT<i>b RN \ p=0 u (4.44) 

The meson-nucleon couplings are related to the 1PI vertices by 

(N\r] a N) = g v <*NN uj 5 u = iuT va fj N u (4.45) 

Re-expressing eq.(4.44) in terms of the canonically normalised 'OZI boson' field r) a , we 
therefore derive 

2m N G a A (0) = f aa g rNN (4.46) 

This relation will be useful to us when we consider the 'proton spin' problem. 

All that now remains to cast this into its final form is to make the familiar change 
of variables from Q,r) a to G,r/ a , where rf 1 are interpreted as the physical mesons. We 
therefore find the generalised U(1)a Goldberger-Treiman relation: 

2m N G a A (0) = f aa g v *NN + ^j2^]Ag GNN 5 a0 (4.47) 

For the individual components, this is 

2m N G\ = f n g nNN (4.48) 
2m N G\ = f 8v 'g v > NN + f v g V NN (4.49) 
2m N G° A = f* grfNN + f 0l 'g v NN + V6Ag GNN (4.50) 

The renormalisation group properties of these relations are described in great detail in 
ref.[9]. It is clear that the flavour singlet axial coupling G A satisfies a homogeneous RGE 
and scales with the anomalous dimension 7 corresponding to the multiplicative renormal- 
isation of J° 5 . In the form (4.46), RG consistency is simply achieved by 

Vf aa = jS a0 f aa Vg rNN = (4.51) 

All the scale dependence is in the decay constant f 0a while the the coupling g^NN of the 
'OZI boson' to the nucleon is RG invariant (in contrast to g^^)- In the final form (4.47) 
involving the physical decay constants, a careful analysis shows that apart from G° A (0) the 
only other non RG-invariant quantity is the gluonic coupling g G NN, which is required to 
satisfy the following non-homogeneous RGE to ensure the self-consistency of eq.(4.50): 

Vg GNN = j(g GN N+ ^— ^f a g v <*NN) (4.52) 



8 The p — » limit is delicate, as is the case for the derivation of the conventional Goldberger-Treiman 
relation, and should be taken in this order. Literally at p = 0, both sides vanish since u-y^u = 0. 
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The large-iV c behaviour in the flavour singlet relation is as follows: G\ = 0(N C ), 
f V j0v> = {y/TQ, A = 0(1), g v NN,gr,>NN = 0{^fW c ), g GNN = O(l). So the final 
term AgcNN is 0(1), suppressed by a power of 1/N C compared to all the others, which are 
0(N C ). 

We see that, like gayy, the gluonic coupling gcNN is suppressed at large N c relative 
to the corresponding meson couplings. However, unlike gcyy which is RG invariant, gcNN 
has a complicated RG non-invariance and depends on the anomaly-induced anomalous 
dimension 7. The conjecture in section 4.3 then suggests that while the OZI or large-iV c 
approximation should be a good guide to the value of 5^775 we may expect significant 
OZI violations for gGNN- We would therefore not be surprised to find that gcNN makes a 
sizeable numerical contribution to the U(1)a Goldberger-Treiman relation. 

We now try to test these expectations against the experimental data. We first introduce 
a notation that has become standard in the literature on deep-inelastic scattering. There, 
the axial couplings are written as 

G\ = - a 3 G\ = —^-= a 8 G° A = a (4.53) 
A 2 A 2^/3 A VQ 

where the a a have a simple interpretation in terms of parton distribution functions. 
Experimentally, 

a 3 = 1.267 ±0.004 a 8 = 0.585 ± 0.025 (4.54) 

from low-energy data on nucleon and hyperon beta decay. The latest result 9 for a quoted 
by the COMPASS collaboration [53] from deep-inelastic scattering data is 

a | Q2 _ +oo = 0.33 ±0.06 (4.55) 

with a similar result from HERMES [54] . 

The OZI expectation is that a = a 8 . In the context of DIS, this is a prediction of the 
simple quark model, where it is known as the Ellis- Jaffe sum rule [57]). We return to this 
in the context of the 'proton spin' problem in section 5 but for now we concentrate on the 
low-energy phenomenology of the pseudoscalar meson-nucleon couplings. 

The original Goldberger-Treiman relation (4.48) gives the following value for the pion- 
nucleon coupling, 

g nNN = 12.86 ±0.06 (4.56) 

consistent to within about 5% with the experimental value 13.65(13.80) ±0.12 (depending 
on the dataset used [58]). In an ideal world where g^NN and g^'NN were both known, 
we would now verify the octet formula (4.49) then determine the gluonic coupling gcNN 
from the singlet Goldberger-Treiman relation (4.50). However, the experimental situation 
with the i] and r/-nucleon couplings is far less clear. One would hope to determine these 



9 This supersedes the result a°| Q 2 =4GeV 2 = 0.2371^ quoted by COMPASS in 2005 [55, 56], which we 
used as input into our analysis of the phenomenology of the U(1)a GT relation in ref.[3]. The fits presented 
here are updated from those of ref.[3] to take account of this. For a further discussion of the experimental 
situation, see section 5. 
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Figure 6: These figures show the dimensionless 77-nucleon coupling g V NN and the gluonic coupling 
9gnn in units of GeV -3 expressed as functions of the experimentally uncertain r/-nucleon coupling 
9ri'NN, as determined from the flavour octet and singlet Goldberger-Treiman relations (4.49) and 
(4.50). 

couplings from the near threshold production of the rj and rf in nucleon-nucleon collisions, 
i.e. pp — > pprj and pp — > pprj' , measured for example at COSY-II [59, 60, 61]. However, 
the r\ production is dominated by the N(1535)Su nucleon resonance which decays to Nrj, 
and as a result very little is known about g v NN itself. The detailed production mechanism 
of the rf is not well understood. However, since there is no known baryonic resonance 
decaying into Nrj', we may simply assume that the reaction pp — * ppr( is driven by the 
direct coupling supplemented by heavy-meson exchange. This allows an upper bound to 
be placed on g v >NN and on this basis ref.[62] quotes g v >NN < 2.5. This is supported by an 
analysis [63] of very recent data from CLAS [64] on the photoproduction reaction "fp — > prf . 
Describing the cross-section data with a model comprising the direct coupling together with 
i-channel meson exchange and s and u-channel resonances, it is found that equally good 
fits can be obtained for several values of g n 'NN covering the whole region < g v 'NN < 2.5. 

In view of this experimental uncertainty, we shall use the octet and singlet Goldberger- 
Treiman relations to plot the predictions for g^NN and gcNN as a function of the ill- 
determined r/-nucleon coupling in the experimentally allowed range < g^NN < 2.5. The 
results (again taking the value (4.32) for A) are given in Fig. 6. In Fig. 7 we have shown 
the relative magnitudes of the various contributions to the flavour-singlet formula. 

What we learn from this is that for values of g^NN approaching the upper end of 
the experimentally allowed range, the contribution of the OZI-suppressed gluonic coupling 
gGNN is quite large. The variation of f 0r> g^'NN over the allowed range is compensated 
almost entirely by the variation of VQgcNN, with the f^g^NN contribution remaining 
relatively constant. 

For example, if experimentally we found gwNN — 2-5, which corresponds to the cross- 
sections for pp — > ppr( and 7p — > prf being almost entirely determined by the direct 
coupling, then we would have g v NN — 4.14 and goNN — —31.2 GeV -3 . In terms of the 
contributions to the U(1)a Goldberger-Treiman relation, this would give (in GeV) 

2m N G° A [N c ; 0.25] = /°"' 9rfNN [N c ; 0.26] + f 0r >g V NN[N c ; 0.09] 

+ V6Ag G NN[0(l); -0.10] (4.57) 

The anomalously small value of G\ compared to its OZI value (the OZI approximation is 
2mArG^| OZI = y/2 2m^G\ = 0.45) is then due to the partial cancellation of the sum of the 
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Figure 7: This shows the relative sizes of the contributions to the U(1)a Goldberger-Treiman 
relation from the individual terms in eq.(4.50), expressed as functions of the coupling g v >NN- The 
dotted (black) line denotes 2itinG a . The long-dashed (green) line is f° n g n >NN and the short- 
dashed (blue) line is f° n g V NN- The solid (red) line shows the contribution of the novel gluonic 
coupling, \/6AgcNN, where A determines the QCD topological susceptibility. 

meson-nucleon coupling terms by the gluonic coupling gcNN- Although formally 0(1/N C ) 
suppressed, numerically it gives a major contribution to the large OZI violation in Ga. 
This would give some support to our conjecture and provide further evidence that we are 
able to predict the location of large OZI violations using the renormalisation group as a 
guide. 

Of course, it may be that experimentally we eventually find a value for g v 'NN — 1-5, 
in the region where gGNN contributes only around 10% or less. Although surprising, this 
would open the possibility that all gluonic couplings of type gcxx are close to zero, which 
could be interpreted as implying that the gluonic component of the rf wave function is 
simply small. Clearly, a reliable determination of g v >NN, or equivalently g v NN, would shed 
considerable light on the U{\)a dynamics of QCD. 

5. Topological charge screening and the 'proton spin' 

So far, we have focused on the implications of the U{\)a anomaly for low-energy QCD 
phenomenology. However, the anomaly also plays a vital role in the interpretation of 
high-energy processes, in particular polarised deep-inelastic scattering. 

In this section, we discuss one of the most intensively studied topics in QCD of the 
last two decades - the famous, but misleadingly named, 'proton spin' problem. We review 
the interpretation initially proposed by Veneziano [4] and developed by us in a series of 
papers exploring the relation with the U(1)a GT relation and gluon topology [8, 9, 65]. 
In subsequent work with Narison, we were able to quantify our prediction by using QCD 
spectral sum rules to compute the slope x'(0) of the topological susceptibility [10, 52]. 
Remarkably, the most recent experimental data from the COMPASS [53] and HERMES 
[54] collaborations, released in September 2006, now confirms our original 1994 numerical 
prediction [10]. 




2 . 5 
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5.1 The and angular momentum sum rules 

The 'proton spin' problem concerns the sum rule for the first moment of the polarised 
proton structure function gf. This is measured in polarised DIS experiments through the 
inclusive processes fip -► fiX (EMC, SMC, COMPASS at CERN) or ep -► eX (SLAC, 
HERMES at DESY) together with similar experiments on a deuteron target. The polari- 
sation asymmetry of the cross-section is expressed as 

dAa Y P WTT 2 a 2 /MV" 



'-9l(x,Q 2 ) + 0(^-) (5-1) 



dxdy 2 « \ Q 2 

with conventional notation: Q 2 = — q 2 and x = Q 2 /2p 2 .q are the Bjorken variables, where 
P2, q are the momenta of the target proton and incident virtual photon respectively, y = 
Q 2 /xs and Y p = (2 - y)/y. 

According to standard theory, g\ is determined by the proton matrix element of two 
electromagnetic currents carrying a large spacelike momentum. The sum rule for the first 
moment of g\ is derived from the twist 2, spin 1 terms in the operator product expansion 
for the currents: 

J\q)n-q) qZoo 2^|i [ACf(« s )(^ 5 + -LjJ,) + H^ACf [a s )J%] (5.2) 

where AC^ S and ACf are Wilson coefficients and J^ 5 (a = 3, 8, 0) are the renormalised 
axial currents, with the normalisations defined in section 2. It is the occurrence of the axial 
currents in this OPE that provides the link between the U(1)a anomaly and polarised DIS. 
The sum rule is therefore: 

T?(Q 2 ) ee j\x g^Q 2 ) = ^ACf^a 3 + \a*) + ±AC?a°(Q 2 ) (5.3) 

where the axial charges a 3 , a 8 and a°(Q 2 ) are defined in terms of the forward proton matrix 
elements as in eq.(4.53). Here, we have explicitly shown the Q 2 scale dependence associated 
with the RG non-invariance of a°(Q 2 ). 

Since the flavour non-singlet axial charges are known from low-energy data, a measure- 
ment of the first moment of g\ amounts to a determination of the flavour singlet a°(Q 2 ). At 
the time of the original EMC experiment in 1988 [66] the theoretical expectation based on 
the quark model was that a = a 8 . The resulting sum rule for g\ is known as the Ellis- Jaffe 
sum rule [57]. The great surprise of the EMC measurement was the discovery that in fact a 
is significantly suppressed relative to a 8 , and indeed the earliest results suggested it could 
even be zero. However, the reason the result sent Shockwaves through both the theoretical 
and experimental communities (to date, the EMC paper has over 1300 citations) was the 
interpretation that this implies that the quarks contribute only a fraction of the total spin 
of the proton. 

In fact, this interpretation relies on the simple valence quark model of the proton and 
is not true in QCD, where the axial charge decouples from the real angular momentum 
sum rule for the proton. Rather, as we shall show, the suppression of a°(Q 2 ) reflects the 
dynamics of gluon topology and appears to be largely independent of the structure of the 
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proton itself. Precisely, it is a manifestation of topological charge screening in the QCD 
vacuum. 



The angular momentum sum rule is derived by taking the forward matrix element of 
the conserved angular momentum current M^ V ' K , defined in terms of the energy- momentum 
tensor as 

M ixv\ = x W T ^ + d pX P^ x (5.4) 

The inclusion of the arbitrary tensor X ppuX just reflects the usual freedom in QFT of 
defining conserved currents. This gives us some flexibility in attempting to write M pvX as 
a sum of local operators, suggesting interpretations of the total angular momentum as a 
sum of 'components' of the proton spin. In fact, however, it is not possible to write M^ vX 
as a sum of operators corresponding to quark and gluon spin and angular momentum in a 
gauge-invariant way. The best decomposition is [67, 68, 69] 

Mi xv\ = q^x + q^ x u] + ^ x v] + _ _ _ (5_ 5 ) 
where the dots denote terms whose forward matrix elements vanish. Here, 

Of " = l^q lal5 q = \e^^2n- f J^ 
0% X = iff D X q 

M =F np F X ( 5 _ 6 ) 

At first sight, Of A looks as if it could be associated with 'quark spin', since for free 
Dirac fermions the spin operator coincides with the axial vector current. O^x^ would 
correspond to 'quark orbital angular momentum', leaving 0^ X x^ as 'gluon total angular 
momentum'. Any further decomposition of the gluon angular momentum is necessarily not 
gauge invariant. 

The forward matrix elements of these operators may be expressed in terms of form 
factors and, as we showed in ref.[68], this exhibits an illuminating cancellation. After some 
analysis, we find: 

(p,s\0^ X \p,s) = a°m N e^s a 

(p, s|0£ [ Vl \p, s) = J q -^—p p p^e^P a s a - a°m N e^ XcT s CT 

(p, S \0fx^\p,s) = J g J— Pp p{^ x ^s a (5.7) 
The angular momentum sum rule for the proton is then just 

\ = J q + J 9 (5-8) 

where the Lorentz and gauge-invariant form factors J q and J g may reasonably be thought 
of as representing quark and gluon total angular momentum. However, even this interpre- 
tation is not at all rigorous, not least because J q and J g mix under renormalisation and 
scale as 



dlnQ 2 \J g J 4vr V |Cf -3"// V J g 
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Only the total angular momentum is Lorentz, gauge and scale invariant. 10 

The crucial observation, however, is that the axial charge aP explicitly cancels from 
the angular momentum sum rule, a is an important form factor, which relates the first 
moment of g\ to gluon topology via the U(1)a anomaly, but it is not part of the angular 
momentum sum rule for the proton. 

Just as a can be measured in polarised inclusive DIS, the form factors J q and J g can 
be extracted from measurements of unpolarised generalised parton distributions (GPDs) in 
processes such as deeply-virtual Compton scattering 7*p — > jp. These can also in principle 
be calculated in lattice QCD. The required identifications with GPDs are given in ref.[68]. 

5.2 QCD parton model 

Before describing our resolution of the 'proton spin' problem, we briefly review the parton 
model interpretation of the first moment sum rule for g^. 

In the simplest form of the parton model, the proton structure at large Q 2 is described 
by parton distributions corresponding to free valence quarks only. The polarised structure 
function is given by 

t=i 

where A<fe(x) is defined as the difference of the distributions of quarks (and antiquarks) 
with helicities parallel and antiparallel to the nucleon spin. It is convenient to work with 
the conventionally-defined flavour non-singlet and singlet combinations Aq NS and Aq s 
(often also written as AS). 

In this model, the first moment of the singlet quark distribution Aq s = dx Aq s (x) 
can be identified as the sum of the helicities of the quarks. Interpreting the structure 
function data in this model then leads to the conclusion that the quarks carry only a 
small fraction of the spin of the proton. There is indeed a real contradiction between the 
experimental data and the free valence quark parton model. 

However, this simple model leaves out many important features of QCD, the most 
important being gluons, RG scale dependence and the chiral C/a(1) anomaly. When these 
effects are included, in the QCD parton model, the naive identification of Aq s with spin 
no longer holds and the experimental results for g\ can be accommodated, though not 
predicted. 

In the QCD parton model, the polarised structure function is written in terms of both 
quark and gluon distributions as follows: 

9 P 1 (x,Q 2 ) = 

/;^I[ac-(3a^ ( m) + a^Qa/(^ ) + a^(3a,(m) 



(5.11) 



10 For a careful discussion of the parton interpretation of longitudinal and transverse angular momentum 
sum rules, see ref.[70]. This confirms our assertion that the axial charge a is not to be identified with 
quark helicities in the parton model. 
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where AC S , AC 9 and AC NS are perturbatively calculable functions related to the Wilson 
coefficients and the quark and gluon distributions have a priori a t = lnQ 2 /A 2 dependence 
determined by the RG evolution, or DGLAP, equations. The first moment sum rule is 
therefore 



1 



T\(Q 2 ) = ^ [ACf *Aq NS + ACf Aq s + ACf A 5 j (5.12) 

Comparing with eq.(5.3), we see that the axial charge a°(Q 2 ) is identified with a linear 
combination of the first moments of the singlet quark and gluon distributions. It is often, 
though not always, the case that the moments of parton distributions can be identified 
in one-to-one correspondence with the matrix elements of local operators. The polarised 
first moments are special in that two parton distributions correspond to the same local 
operator. 

The RG evolution equations for the first moments of the parton distributions are 
derived from the matrix of anomalous dimensions for the lowest spin, twist 2 operators. 
This introduces an inevitable renormalisation scheme ambiguity in the definitions of Aq 
and Ag, and their physical interpretation is correspondingly nuanced. The choice closest to 
our own analysis is the 'AB' scheme [71] where the parton distributions have the folowing 
RG evolution: 

u^^ NS = o z4^ s = 

j£c?tA9(Q 2 ) = l{^9(Q 2 ) - |Ag*) (5.13) 

which requires ACf = . It is then possible to make the following identifications 

with the axial charges: 

a 3 = An - Ad 

a 8 = Au + Ad - 2As 

3(T 

a°(Q 2 ) = Au + Ad + As- — -Ag(Q 2 ) (5.14) 

with Aq s = Au + Ad+ As. Notice that in the AB scheme, all the scale dependence of the 
axial charge a°(Q 2 ) is assigned to the gluon distribution Ag(Q 2 ). 

This was the identification originally introduced for the first moments by Altarelli and 
Ross [72], and resolves the 'proton spin' problem in the context of the QCD parton model. 
In this scheme, the Ellis-Jaffe sum rule follows from the assumption that in the proton both 
As and Ag(Q 2 ) are zero, which is the natural assumption in the free valence quark model. 
This is equivalent to the OZI approximation a°(Q 2 ) = a 8 . However, in the full QCD parton 
model, there is no reason why Ag(Q 2 ), or even As, should be zero in the proton. Indeed, 
given the different scale dependence of a°(Q 2 ) and a 8 , it would be unnatural to expect this 
to hold in QCD itself. 

An interesting conjecture [72] is that the observed suppression in a°(Q 2 ) is due over- 
whelmingly to the gluon distribution Ag(Q 2 ) alone. Although by no means a necessary 
consequence of QCD, this is a reasonable expectation given that it is the anomaly (which 
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is due to the gluons and is responsible for OZI violations) which is responsible for the scale 
dependence in a°(Q 2 ) and Ag(Q 2 ) whereas the Aq are scale invariant. This would be in 
the same spirit as our conjecture on OZI violations in low-energy phenomenology in section 
4.3. To test this, however, we need to find a way to measure Ag(Q 2 ) itself, rather than the 
combination a°(Q 2 ). The most direct option is to extract Ag(x,Q 2 ) from processes such 
as open charm production, 7*3 — > cc, which is currently being intensively studied by the 
COMPASS [73], STAR [74] and PHENIX [75] collaborations. 

5.3 Topological charge screening 

We now describe a less conventional approach to deep-inelastic scattering based entirely 
on field-theoretic concepts. In particular, the role of parton distributions is taken over by 
the 1PI vertices of composite operators introduced above. (For a review, see ref.[76]). 

Once again, the starting point is the use of the OPE to express the moments of a 
generic structure function F(x,Q 2 ) as 

f dxx n ~ l F(x,Q 2 ) = ^2°a(Q 2 ) (p\O n A (0)\p) (5.15) 
Jo A 

where 0\ denotes the set of lowest twist, spin n operators and C\(Q 2 ) are the correspond- 
ing Wilson coefficients. The next step is to introduce a new set of composite operators 
Ob, chosen to encompass the physically relevant degrees of freedom, and write the matrix 
element as a product of two-point Green functions and 1PI vertices as follows: 

1 dx x^Fix, Q 2 ) = C a(Q 2 ) (0\T 0\ 6 b \0) T 6Bpp (5.16) 

A B 

This decomposition splits the structure function into three parts - first, the Wilson coef- 
ficients C\(Q 2 ) which can be calculated in perturbative QCD; second, non-perturbative 
but target independent Green functions which encode the dynamics of the QCD vacuum; 
third, non-perturbative vertex functions which characterise the target by its couplings to 
the chosen operators Ob- 11 

Now specialise to the first moment sum rule for g\. For simplicity, we first present the 
analysis for the chiral limit, where there is no flavour mixing. Using the anomaly (2.4), we 
can express the flavour singlet contribution to the sum rule as 

F?(Q 2 ) S i ng i e t = J dx g\{x,Q 2 ) singlet = ACf (a s ) ( P \Q\ P ) (5.17) 

The obvious choice for the operators ($b in this case are the flavour singlet pseudoscalars 
and it is natural to choose the 'OZI boson' field fp = f 00 -^^, which is normalised so 

11 Wc emphasise again that this decomposition of the matrix elements into products of Green functions 
and 1PI vertices is exact, independent of the choice of the set of operators Ob- In particular, it is not 
necessary for Ob to be in any sense a complete set. If a different choice is made, the vertices Tq bPP 
themselves change, becoming 1PI with respect to a different set of composite fields. In practice, the set 
of operators Ob should be as small as possible while still capturing the essential degrees of freedom. A 
good choice can also result in vertices ^o BVV which are both RG invariant and closely related to low-energy 
physical couplings. 
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Figure 8: Illustration of the decomposition of the matrix element (p\Q\p) into two-point Green 
functions and 1PI vertices. The Green function in the first diagram is x(0); m the second it is 



that d/dp 2 r^o| = 1. As we have seen in eq.(4.51), the corresponding 1PI vertex is 
then RG invariant. Writing the 1PI vertices in terms of nucleon couplings as in eq.(4.45), 
we find 



r?(Q 2 : 



singlet 



H^i-ACf K) (<0|T Q Q\0) g QNN + (0\T Q fj°\0) ^ 0AW ) (5.18) 



Recalling that the matrix of two-point Green functions is given by the inversion formula 



Wee Wes. 
Ws, e Ws^ s, 



T QQ F Qr)° 



(5.19) 



and using the normalisation condition for fj°, we can easily show that at zero momentum, 



it; 



es. 



^W ee \ p=0 



(5.20) 



Finally, therefore, we can represent the first moment of g\ in the following, physically 
intuitive form: 



I | (Q ) singlet 



2 1 



3 2mjv 



ACf{a s ) ( X (Q) gQNN + 9f,o 



NN 



(5.21) 



This shows that the first moment is determined by the gluon topological susceptibility in 
the QCD vacuum as well as the couplings of the proton to the pseudoscalar operators Q 
and fj°. In the chiral limit, x(0) = so the first term vanishes. The entire flavour singlet 
contribution is therefore simply 



) singlet 



2 1 



3 2m n 



ACf (a s ) v9(0) 9fp NN 



(5.22) 



The 1PI vertex g^o^N ls invariant, and we see from eq.(2.25) that in the chiral limit 
the slope of the topological susceptibility scales with the anomalous dimension 7, viz. 



d 



(5.23) 



cilnQ 2 

ensuring consistency with the RGE for the flavour singlet axial charge. 

The formulae (5.21) and (5.22) are our key result. They show how the first moment 
of g\ can be factorised into couplings gQNN and g^NN which carry information on the 
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proton structure, and Green functions which characterise the QCD vacuum. In the case of 
eft, the Green functions reduce simply to the topological susceptibility x(0) and its slope 
x'(0). We now argue that the experimentally observed suppression in the first moment 
of gl is due not to a suppression in the couplings, but to the vanishing of the topological 
susceptibility x(0) and an anomalously small value for its slope x'(0). This is what we 
refer to as topological charge screening in the QCD vacuum. 

The justification follows our now familiar conjecture on the relation between OZI vio- 
lations and RG scale dependence. We expect the source of OZI violations to be in those 
quantities which are sensitive to the anomaly, as identified by their scaling dependence on 
the anomalous dimension 7, in this case x'(0). In contrast, it should be a good approxima- 
tion to use the OZI value for the RG- invariant vertex ffjjOjvjV) that is gffiNN — V^9fi»NN- 
The corresponding OZI value for \J x'(0) would be / 7r /v / 6. This gives our key formula for 
the flavour singlet axial charge: 

The corresponding prediction for the first moment of eft is 

T\(Q 2 ) smg iet = \^Cf{a s ) a 8 ^ y^W) (5.25) 

The final step is to compute the slope of the topological susceptibility. In time, lattice 
gauge theory should provide an accurate measurement of x'(0)- However, this is a par- 
ticlarly difficult correlator for lattice methods since it requires a simulation of QCD with 
light dynamical fermions and algorithms that implement topologically non-trivial configu- 
rations in a sufficiently fast and stable way. Instead, we have estimated the value of x'(0) 
using the QCD spectral sum rule method. Full details and discussion of this computation 
can be found in refs.[10, 52]. The result is: 

Vx'(O) = 26.4 ± 4.1 MeV (5.26) 

This gives our final prediction for the flavour singlet axial charge and the complete 
first moment of eft: 

aC V=iOGev2 = 0.33 ±0.05 (5.27) 
r ?lcP=iOGev2 = 0.144 ±0.009 (5.28) 

Topological charge screening therefore gives a suppression factor of approximately 0.56 in 
a compared to its OZI value a 8 = 0.585. 

In the decade since we made this prediction, the experimental measurement has been 
somewhat lower than this value, in the range a ~ 0.20 — 0.25. This would have suggested 
there is also a significant OZI violation in the nucleon coupling g^o^N itself, implicating the 
proton structure in the anomalous suppression of P^. Very recently, however, the COM- 
PASS and HERMES collaborations have published new results on the deuteron structure 
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Figure 9: COMPASS and SMC data for the deuteron structure function gf(x). Statistical error 
bars are shown with the data points. The shaded band shows the systematic error. 



function which spectacularly confirm our picture that topological charge screening in the 
QCD vacuum is the dominant suppression mechanism. 

This new data is shown in Fig. 9. This is based on data collected by COMPASS 
at CERN in the years 2002-2004 and has only recently been published. The accuracy 
compared to earlier SMC data at small x is significantly improved and the dip in xgf 
around x ~ 10~ 2 suggested by the SMC data is no longer present (Fig. 9). This explains 
the significantly higher value for a found by COMPASS compared to SMC. From this 
data, COMPASS quote the first moment for the proton-neutron average g^ = {g\ + <7™)/2 
as [53] 

rf | Q 2 =3GeV 2 = 0.050 ± 0.003(stat) ± 0.002(euoZ) ± 0.005(syst) (5.29) 
Extracting the flavour singlet axial charge from the analogue of eq.(5.3) for then gives 

a°| Q2=3GeV 2 = 0.35 ± Om(stat) ± 0.05{syst) (5.30) 

or evolving to the Q 2 — > oo limit, 

a | Q 2_ oo = 0.33 ± 0.03(stat) ±0.05(syst) (5.31) 

Similar results are found by HERMES, who quote [54] 

a°| Q 2 =5GeV 2 = 0.330 ± 0.011(th) ± 0.025(exp) ±0.028(evol) (5.32) 

The agreement with our prediction (5.27) is striking. 

To close this section, we briefly comment on the extension of our analysis beyond the 
chiral limit. In this case, the operator ^/2njQ in eq.(5.17) is replaced by the full divergence 
of the flavour singlet axial current, viz. D° = sJIn^Q + dob c m b <p^. Separating the matrix 
element (p\D°\p) into Green functions and 1PI vertices, we find from the zero-momentum 
Ward identities that (0\T D° Q\0) = so the contribution from gQNN still vanishes. The 
other Green function is (0\T D° fj a \0) = —f 0a , so the first moment sum rule becomes 

T^Q 2 ) singlet = J r^- ACf (a s ) f^g^ NN (5.33) 
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It is clear that this is simply an alternative derivation of the £7(1) GT relation (4.46) for 
a . We could equally use the alternative form (4.47) to write 

r p 1 (Q 2 ) Sin9 iet = AC?(a s ) v/6 (7°%™ + \fcAg GNN ) (5.34) 

Recalling the RGE (4.52) for gcNN, we see that this bears a remarkable similarity to the 
expression for a in terms of parton distributions in the AB scheme, eq.(5.14). This was 
first pointed out in ref.[8, 9]. 

Manipulating the zero-momentum Ward identities in a similar way to that explained 
above in the chiral limit now shows that we can express the decay constants f aa in terms 
of vacuum Green functions as follows (see eq.(3.9): 

UFU = ^(0\T D a D b \Q)\ p=0 (5.35) 

However, for non-zero quark masses there is flavour mixing amongst the 'OZI bosons' rj a 
and we cannot extract the decay constants simply by taking a square root, as was the 
case in writing / 00 = \/x'(fy in the chiral limit. Nevertheless, in ref.[52] we estimated the 
decay constants and form factors in the approximation where we use eq.(5.35) with the full 
divergence D a but neglect flavour mixing. Assuming OZI for the couplings, this gives the 
estimate 

where we take 



?oo 



(5.36) 



/ 00 ^ ^~(0\T D° D°\0)\ p=0 /» ~ ^~{0\T D* D*\0)\ p=0 (5.37) 

Evaluating the Green functions using QCD spectral sum rules gives 

a( V=iOGcv2 = 0.31 ±0.02 (5.38) 
r ilQ 2 =iOGcV^ = 0.141 ±0.005 (5.39) 

As we have seen in the last section, flavour mixing can be non-negligible in the phenomenol- 
ogy of the pseudoscalar mesons, so we should be a little cautious in over-estimating the 
accuracy of these estimates. (The quoted errors do not include this systematic effect.) 
Nevertheless, the fact that they are consistent with those obtained in the chiral limit re- 
inforces our confidence that the flavour singlet axial charge is relatively insensitive to the 
quark masses and that eqs.(5.27) and (5.28) indeed provide an accurate estimate of the 
first moment of g\. 

The observation that the 'proton spin' sum rule could be explained in terms of an 
extension of the Goldberger-Treiman relation to the flavour singlet sector was made in 
Veneziano's original paper [4]. This pointed out for the first time that the suppression 
in a was an OZI-breaking effect. Since the Goldberger-Treiman relation connects the 
pseudovector form factors with the pseudoscalar channel, where it is known that there are 
large OZI violations for the flavour singlet, it becomes natural to expect similar large OZI 
violations also in a . This is the fundamental intuition which we have developed into a 
quantitative resolution of the 'proton spin' problem. 
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Figure 10: Semi-inclusive DIS eN — > ehX in the target fragmentation region. In the equivalent 
current fragmentation process, the detected hadron h is emitted from the hard collision with 7. 
The right-hand figure shows a simple Reggeon exchange model valid for z ~ 1, where h carries a 
large target energy fraction. 

5.4 Semi-inclusive polarised DIS 

While the agreement between our prediction for the first moment of g\ and experiment is 
now impressive, it would still be interesting to find other experimental tests of topological 
charge screening. A key consequence of this mechanism is that the OZI violation observed 
in a is not a property specifically of the proton, but is target independent. This leads us to 
look for ways to make measurements of the polarised structure functions of other hadronic 
targets besides the proton and neutron. We now show how this can effectively be done by 
studying semi- inclusive DIS eN — > ehX in the target fragmentation region (Fig. 10). 

The differential cross-section in the target fragmentation region can be written analo- 
gously to eq.(5.1) in terms of fracture functions: 



where x = Q 2 /2p2.q, xb = Q 2 /2k.q, z = P2-q/p2-q so that 1 — z = x/xb, and the invariant 
momentum transfer t = K 2 = — k 2 , where k is the momentum of the struck parton. For 
K 2 <C Q 2 , z — Eh/ En (in the photon- nucleon CM frame) is the energy fraction of the 
target nucleon carried by the detected hadron h. 

AMi N is the fracture function [77] equivalent of the inclusive structure function g± , 
so in the same way as in eq.(5.10) we have 



Here, AMj lN (x, z,t, Q 2 ) is an extended fracture function, introduced by Grazzini, Trentadue 
and Veneziano [78], which carries an explicit dependence on t. One of the advantages of 
these fracture functions is that they satisfy a simple, homogeneous RG evolution equation 
analogous to the usual inclusive parton distributions. 

Our proposal [11, 79] (see also [80]) is to study semi-inclusive DIS in the kinematical 
region where the detected hadron h (tt, K or D) carries a large target energy fraction, i.e. 
z approaching 1, with a small invariant momentum transfer t. In this region, it is useful 




(5.40) 




(5.41) 
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to think of the target fragmentation process as being simply modelled by a single Reggeon 
exchange (see Fig. 10), i.e. 

AM? N (x,z,t,Q 2 )\ z ^ ~ F(t)(l-z)- 2a ^gf(x B ,t,Q 2 ) (5.42) 

If we consider ratios of cross-sections, the dynamical Reggeon emission factor F(t)(l — 
z ^-2ot B (t) cance i anc i we -vvill be able to isolate the ratios of gf(xts,t,Q 2 ) for different 
effective targets B. Although single Reggeon exchange is of course only an approximation 
to the more fundamental QCD description in terms of fracture functions (see ref.[81] for 
a more technical discussion), it shows particularly clearly how observing semi-inclusive 
processes at large z with particular choices of h and N amounts in effect to performing 
inclusive DIS on virtual hadronic targets B. Since our predictions will depend only on the 
SU (3) properties of B, together with target independence, they will hold equally well when 
B is interpreted as a Reggeon rather than a pure hadron state. 

The idea is therefore to make predictions for the ratios 1Z of the first moments of the 
polarised fracture functions Jq z dxAM± N (x, z, t, Q 2 ) or equivalently Jq 1 dx^gf(xts,t,Q 2 ) 
for various reactions. The first moments Ff are calculated as in eq.(5.3) in terms of the 
axial charges a 3 , a 8 and a°(Q 2 ) for a state with the SU(3) quantum numbers of B. We then 
use topological charge screening to say that a°(Q 2 ) ~ s(Q 2 )a°| OZI , i.e. the flavour singlet 
axial charge is suppressed relative to its OZI value by a universal, target-independent, 
suppression factor s(Q 2 ). From our calculation of \J %' (0) and the experimental results for 
<7i, we have s|q2 =10 g c v 2 — 0.33/0.585 = 0.56. 

Some of the more interesting predictions obtained in ref.fll] are as follows. The ratio 

/en -► evr+AX 2s - 1 

K[ ~ (5.43) 

\ep-Knr-X J gnjl 2s + 2 y ' 

is calculated by comparing Y\ for the A - and A ++ . It is particularly striking because 
the physical value of s(Q 2 ) is close to one half, so the ratio becomes very small. For 
strange mesons, on the other hand, the ratio depends on whether the exchanged object 
is in the 8 (where the reduced matrix elements involve the appropriate F/D ratio) or 10 
representation, so the prediction is less conclusive, viz. 

„(en-*eK+X\ 2s - 1 - 3(2s - \)F/D 2s - 1 , . . „ 

* U-« " 2.-1 -3(2. +i)f/d (8) or 57TT (10) (5 - 44) 

which we find by comparing Ti for either the T," and S + in the 8 representation or X*~ 
and S* + in the 10. For charmed mesons, we again find 

^fen^eD°X\ 2s - 1 

n[ — - ~ (5.45) 

\ep^eD-Xj z ^ 2s + 2 v ; 

corresponding to the ratio for S° to 

At the other extreme, for z approaching 0, the detected hadron carries only a small 
fraction of the target nucleon energy. In this limit, the ratio TZ of the fracture function 
moments becomes simply the ratio of the structure function moments for n and p, i.e. using 
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the current experimental values, 1Z z ~o ~ T^/T^ = —0.30. This is to be compared with 
the corresponding OZI or Ellis- Jaffe value of —0.12. 

The differences between the OZI, or valence quark model, expectations and our pre- 
dictions based on topological charge screening can therefore be quite dramatic and should 
give a very clear experimental signal. In ref.[79], together with De Florian, we analysed 
the potential for realising these experiments in some detail. Since we require particle iden- 
tification in the target fragmentation region, fixed-target experiments such as COMPASS 
or HERMES are not appropriate. The preferred option is a polarised ep collider. 

The first requirement is to measure particles at extremely small angles (6 < 1 mrad), 
corresponding to t less than around 1 GeV 2 . This has already been achieved at HERA in 
measurements of diffractive and leading proton/neutron scattering using a forward detec- 
tion system known as the Leading Proton Spectrometer (LPS). The technique for measuring 
charged particles involves placing detectors commonly known as 'Roman Pots' inside the 
beam pipe itself. 

The next point is to notice that the considerations above apply equally to p as to 
7r production, since the ratios 1Z are determined by flavour quantum numbers alone. The 
particle identification requirements will therefore be less stringent, especially as the produc- 
tion of leading strange mesons from protons or neutrons is strongly suppressed. However, 
we require the forward detectors to have good acceptance for both positive and negatively 
charged mesons M = it, p in order to measure the ratio (5.43). 

The reactions with a neutron target can be measured if the polarised proton beam 
is replaced by polarised 3 He. In this case, if we assume that 3 He = Ap + Bn, the cross 
section for the production of positive hadrons h + measured in the LPS is given by 

a( 3 He -»• h + ) ~ Aa(p -► h + ) + Ba(n -> p) + Ba(n -► M+) (5.46) 

The first contribution can be obtained from measurements with the proton beam. However, 
to subtract the second one, the detectors must have sufficient particle identification at least 
to distinguish protons from positively charged mesons. 

Finally, estimates of the total rates [79] suggest that around 1% of the total DIS events 
will contain a leading meson in the target fragmentation region where a LPS would have 
non- vanishing acceptance (z > 0.6) and in the dominant domain x < 0.1. The relevant 
cross-sections are therefore sufficient to allow the ratios TZ to be measured. 

The conclusion is that while our proposals undoubtedly pose a challenge to experi- 
mentalists, they are nevertheless possible. Given the theoretical importance of the 'proton 
spin' problem and the topological charge screening mechanism, there is therefore strong 
motivation to perform target fragmentation experiments at a future polarised ep collider 
[82]. 

6. Polarised two-photon physics and a sum rule for g\ 

The U{1)a anomaly plays a vital role in another sum rule arising in polarised deep-inelastic 
scattering, this time for the polarised photon structure function gj(x~,, Q 2 ; K 2 ). For real 
photons, the first moment of gj vanishes as a consequence of electromagnetic current 
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Figure 11: Kinematics for the two-photon DIS process e + e 



e+erX. 



conservation [83]. For off-shell photons, we proposed a sum rule in 1992 [5, 6] whose 
dependence on the virtual momentum of the target photon encodes a wealth of information 
about the anomaly, chiral symmetry breaking and gluon dynamics in QCD. This is of special 
current interest since, given the ultra-high luminosity of proposed e + e~ colliders designed 
as B factories, a detailed measurement of our sum rule is about to become possible for the 
first time. 

6.1 The first moment sum rule for gj 

The polarised structure function gj is measured in the process e + e~ — ► e + e~ X, which at 
sufficiently high energy is dominated by the two-photon interaction shown in Fig. 11. The 
deep-inelastic limit is characterised by Q 2 — > 00 with x = Q 2 /2p2-q and x 7 = Q 2 /2k.q 
fixed, where Q 2 = —q 2 , K 2 = —k 2 and s = (p± +P2) 2 ■ The target photon is assumed to be 
relatively soft, K 2 <C Q 2 . 

We are interested in the dependence of the photon structure function gJ(x 7 ,Q 2 ; K 2 ) 
on the invariant momentum K 2 of the target photon. Experimentally, this is given by 
K 2 ~ EE' 2 6 2 where E' 2 and 9 2 are the energy and scattering angle of the target electron. 
For the values K 2 ~ m 2 of interest in the sum rule, the target electron is nearly-forward 
and 62 is very small. If it can be tagged, then the virtuality K 2 is simply determined 
from 62; otherwise K 2 can be inferred indirectly from a measurement of the total hadronic 
energy. 

The total cross-section a and the spin asymmetry Act can be expressed formally in 
terms of 'electron structure functions' as follows [5] 



where a = 2~(<J++ + cr H ) and Act = 2( <T ++ ~ a -\ ) with +, — referring to the electron 

helicities. The parameter y = Q 2 /xs <C 1 and only the leading order terms are retained 
below. 




(6.1) 




(6.2) 
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These electron structure functions can be expressed as convolutions of the photon 
structure functions with appropriate splitting functions. In particular, we have 

gl(x,Q 2 ) = f r% C^pJ-) 9 ^,Q 2 ;K 2 ) (6.3) 

where AP 7e (x) = (2 — x). This allows us to relate the x 7 -moments of the photon structure 
functions to the x- moments of the cross-sections. For the first moment of gj, we find: 

I dxx dQHx°dK* = r 3 i^l d ^ 9 ^ Q ^ K ^ (6 - 4) 

The first moment sum rule follows, as for the proton, by using the OPE (5.2) to 
express the product of electromagnetic currents for the incident photon in terms of the 
axial currents J^ 5 . The matrix elements {j* (k)\J® 5 \j* (k)) with the target photon are 
then expressed in terms of the 3-current AVV Green function involving one axial and two 
electromagnetic currents. We define form factors for this fundamental correlator as follows: 

-i{0\J^(p)Jx(ki)J p (k 2 )\Q) = A\ €^\p a ki + A% e^ Xpa k% 

+ ^3 e /iAa/3^1 ^2 ^2p + A\ (-ppafik^k^kw 

+ A\ enXapkik^kip + A§ (. jxpoi ^k < ^ k% k 2 \ 

(6.5) 

where the six form factors are functions of the invariant momenta, i.e. Af = Af(p 2 , k 2 ,k%). 
We also abbreviate Af(0, k 2 , k 2 ) = Af(K 2 ). 

The first moment sum rule for gj is then [5]: 

f d Xl gUx^Q 2 ^ 2 ) = Aira £ AC^Q 2 ) (a^K 2 ) - A a 2 {K 2 )) (6.6) 

•* a=3,8,0 

where the Wilson coefficients are related to those in eq.(5.3) by ACf = AC^ 18 , ACf = 
^ACf s and AC? = ^ACf . 12 

Now, just as the sum rule for the proton structure function g\ could be related to 
low-energy meson-nucleon couplings via the U(1)a Goldberger-Treiman relations, we can 
relate this sum rule for gj to the pseudoscalar meson radiative decays using the analysis in 
section 4.2. Introducing the off-shell radiative pseudoscalar couplings for photon virtuality 
K 2 , we define form factors 

F a (K 2 ) = l-^a^y 1 r 9ril (K 2 ) (6.7) 

12 Explicitly, 

at leading order, where t(Q) = \ In Sj- and 7 = — \j^yz is the anomalous dimension corresponding to the 
U(1)a current renormalisation. 
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(6.8) 



where the cig m are the electromagnetic C/(1)a anomaly coefficients defined earlier. We may 
then rewrite the sum rule as 

f dx^gJ(x 7 ,Q 2 ;K 2 ) = \- ^C!(Q 2 )a a cm F a {K 2 ) (6.9) 

J ° Z7T a=3,8,0 

The dependence of the g\ on the invariant momentum K 2 of the target photon reflects 
many key aspects of both perturbative and non-perturbative QCD dynamics. For on-shell 
photons, K 2 = 0, we have simply [83, 5] 

-l 

dx 7 gJ(x 7 ,Q 2 ;K 2 = 0) = (6.10) 



/ 

J o 



This is a consequence of electromagnetic current conservation. This follows simply by 
taking the divergence of eq.(6.5) and observing that in the limit p — ► 0, both A\ and Ai 
are of 0(K 2 ). 13 

In the asymptotic limit where K 2 <C m 2 , a relatively straightforward renormalisation 
group analysis combined with the anomaly equation shows that, for the flavour non-singlets, 
the A\ tend to the value ^fag m . while in the flavour singlet sector, A\ has an additional 
factor depending on the anomalous dimension 7. Using the explicit expressions for the 
Wilson coefficients, we find 
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dx 7 gj (x y , Q 2 ;K 2 <^ m 2 ) 
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(6.11) 



The asymptotic limit is therefore determined by the electromagnetic U(1)a anomaly, with 
logarithmic corrections reflecting the anomalous dimension of the flavour singlet current 
due to the colour U(1)a anomaly. (See also ref.[84] for a NNLO analysis.) 



13 Electromagnetic current conservation in eq.(6.5) implies 

A a 2 = 



Al = A a 3 kl + Al\{p 2 -kl-kl), 



A\k\ + Agio 2 - k\ - k\) 



The chiral limit is special since the form factors can have massless poles and is considered in detail in ref.[6]. 
The sum rule (6.10) still holds. 
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In between these limits, the first moment of gj provides a measure of the form factors 
defining the 3-current AVV Green function, which encodes a great deal of information 
about the dynamics of QCD, especially the non-perturbative realisation of chiral symmetry 
[6]. Equivalently, in the form (6.9), it measures the momentum dependence of the off-shell 
radiative couplings of the pseudoscalar mesons as the form factors F a {K 2 ) vary from to 
1. 

Just as for g^, we can again isolate a dependence on the topological susceptibility 
through the identification of the flavour singlet decay constant / 00 in eq.(6.7) with x'(0) 
in the chiral limit. This time, however, it is unlikely to be a good approximation to set 
the corresponding coupling ^o 77 equal to its OZI value since it is not RG invariant. A 
more promising approximation is to recall from section 4 that the RG invariant gluonic 
coupling gcyy(0) is OZI suppressed and likely to be small. This was confirmed by the 
phenomenological analysis. If we assume this is also true of the off-shell coupling, then 
we may approximate the sum rule for gj entirely in terms of the off-shell couplings of the 
physical mesons 7r°, rj and rj . 

In general, the momentum dependence of the form factors (A\ — A?,) or F a will depend 
on the fermions contributing to the AVV Green function [6] . In the case of leptons, or heavy 
quarks, the crossover scale as the form factors F a (K 2 ) rise from to 1 with increasing K 2 
will be given by the fermion mass. For the light quarks, however, we expect the crossover 
scale to be a typical hadronic scale ~ m p rather than m u ^ s - This can be justified by a 
rough OPE argument and is consistent with old ideas of vector meson dominance [6, 85]. 
This behaviour would be an interesting manifestation of the spontaneous breaking of chiral 
symmetry. 

Once again, therefore, we see a close relation between the realisation of sum rules in 
high-energy deep-inelastic scattering and low-energy meson physics. All these issues are 
discussed at some length in our earlier papers, but here we now turn our attention to the 
vital question of whether the gj sum rule can be measured in current or future collider 
experiments [7]. 

6.2 Cross-sections and spin asymmetries at polarised B factories 

The spin-dependent cross-sections for the two-photon DIS process e + e~ — > e + e~X were 
analysed in refs.[5, 7] taking account of the experimental cuts on the various kinematical 
parameters. Keeping the lower cut on Q 2 as a free parameter, we found the following 
results for the total cross-section and spin asymmetry: 



In order to measure the gj sum rule, we need to find collider parameters such that the spin 
asymmetry is significant in a kinematic region where the total cross-section is still large. 




(6.12) 



and 




(6.13) 
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Figure 12: The left-hand graph shows the total cross-section a (in pb) at SuperKEKB as the 
experimental cut * s varied from 1 to 10 GeV 2 . The right-hand graph shows the spin asymmetry 
Aa/a over the same range of Q^ in . 

A useful statistical measure of the significance of the asymmetry is that \/ La Aa/a S> 1, 
where L is the luminosity. 

When we first proposed the first moment sum rule for gj , the luminosity available from 
the then current accelerators was inadequate to allow it to be studied. For example, for a 
polarised version of LEP operating at s = 10 4 GeV 2 with an annual integrated luminosity 
of L = 100 pb -1 , and optimising the cut at Q min = 10 GeV 2 , we only have a ~ 35 pb and 
Aa/a ~ 0.01. The corresponding annual event rate would be 3.5 x 10 3 and the statistical 
significance \/ La Aa/a ~ 0.5, so even a reliable measurement of the spin asymmetry could 
not be made. 

Clearly, a hugely increased luminosity is required and this has now become available 
with proposals for machines with projected annual integrated luminosities measured in 
inverse attobarns. However, as noted in ref.[5], if this increased luminosity is associated 
with increased CM energy, then the 1/s factor in the spin asymmetry (6.13) sharply reduces 
the possibility of extracting gj. There is also a competition as Q^in * s varied between 
increasing spin asymmetry and decreasing total cross-section. This is particularly evident 
when we analyse the potential of the ILC [86, 87] for measuring the sum rule [7]. We find 
that even optimising the cut, the spin asymmetry is still only of order Act /ct ~ 0.002 
when a itself has fallen to around 15 pb. While, given the high luminosity, this would 
allow a measurement of the first moment of gj integrated over K 2 , a detailed study of the 
-ftT 2 -dependence of the sum rule requires a much greater spin asymmetry. 

This leads us to consider instead the new generation of ultra-high luminosity e + e~ col- 
liders. Although these are envisaged as B factories, these colliders operating with polarised 
beams would, as we now show, be extremely valuable for studying polarisation phenomena 
in QCD. As an example of this class, we take the proposed SuperKEKB collider. (The 
analysis for PEPII is very similar, the main difference being the additional ten-fold increase 
in luminosity in the current SuperKEKB proposals.) 

SuperKEKB is an asymmetric e + e~ collider with s = 132 GeV 2 , corresponding to 
electron and positron beams of 8 and 3.5 GeV respectively. The design luminosity is 
5 x 10 35 cm~ 2 s -1 , which gives an annual integrated luminosity of 5 ab _1 [88]. To see the 
effects of the experimental cut on Q^j n in this case, we have plotted the total cross-section 
and the spin asymmetry in Fig. 12, in the range of Q 2 ^ from 1 to 10 GeV 2 . In this range a 
is falling like 1 / while Aa/a rises to what is actually a maximum at Q^ in = 10 GeV 2 . 
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Taking Q 2 nin = 5 GeV 2 , we find a ~ 12.5 pb with spin asymmetry Aa/a ~ 0.1. The 
annual event rate is therefore 6.25 x 10 7 , with V La Aa/a ~ 750. This combination of a 
very high event rate and the large 10% spin asymmetry means that SuperKEKB has the 
potential not only to measure Act but to access the full first moment sum rule for g[ itself. 
Recall from eq.(6.4) that to measure dx gj(x, Q 2 ; K 2 ) we need not just Aa but the fully 
differential cross-section w.r.t. K 2 as well as x and Q 2 if the interesting non-perturbative 
QCD physics is to be accessed. To measure this, we need to divide the data into sufficiently 
fine K 2 bins in order to plot the explicit K 2 dependence of gj, while still maintaining the 
statistical significance of the asymmetry. The ultra-high luminosity of SuperKEKB ensures 
that the event rate is sufficient, while its moderate CM energy means that the crucial spin 
asymmetry is not overly suppressed by its 1/s dependence. 

Our conclusion is that the new generation of ultra-high luminosity, moderate energy 
e + e _ colliders, currently conceived as B factories, could also be uniquely sensitive to im- 
portant QCD physics if run with polarised beams. In particular, they appear to be the 
only accelerators capable of accessing the full physics content of the sum rule for the first 
moment of the polarised structure function gJ(x,Q 2 ; K 2 ). The richness of this physics, 
in particular the realisation of chiral symmetry breaking, the manifestations of the axial 
U(1)a anomaly and the role of gluon topology, provides a strong motivation for giving 
serious consideration to an attempt to measure the gj sum rule at these new colliders. 
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